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Abstract 

The gauge invariant degrees of freedom of matrix models based on an N x N complex 
matrix, with U(N) gauge symmetry, contain hidden free particle structures. These are 
exhibited using triangular matrix variables via the Schur decomposition. The Brauer 
algebra basis for complex matrix models developed earlier is useful in projecting to a 
sector which matches the state counting of N free fermions on a circle. The Brauer algebra 
projection is characterized by the vanishing of a scale invariant laplacian constructed from 
the complex matrix. The special case of N = 2 is studied in detail: the ring of gauge 
invariant functions as well as a ring of scale and gauge invariant differential operators 
are characterized completely. The orthonormal basis of wavefunctions in this special case 
is completely characterized by a set of five commuting Hamiltonians, which display free 
particle structures. Applications to the reduced matrix quantum mechanics coming from 
radial quantization in M = 4 SYM are described. We propose that the string dual of the 
complex matrix harmonic oscillator quantum mechanics has an interpretation in terms of 
strings and branes in 2 + 1 dimensions. 
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1 Introduction 

1.1 Background and motivations 

There is a class of Gaussian matrix models in D spacetime dimensions 1 x M where the 
two-point function of the matrix up to a trivial spacetime dependence, is 

(4 £t*> = a***. (1.1) 

and where there is an adjoint U(N) action 

Z^gZg\ geU(N). (1.2) 

Our main interest is in D = 4, M = 4 superconformal Yang-Mills (SYM) with gauge 
group U(N) at zero coupling and its dimensional reduction on S 3 to a D = 1 matrix 
harmonic oscillator quantum mechanics. Some of our results apply more generally and in 
particular we make connections to the D = Gaussian complex matrix model considered 
by Ginibre [1]. 

M = 4 super Yang-Mills contains three complex scalar fields in the adjoint of the 
gauge group. The states built from holomorphic functions in one of these fields, say Z, 
comprise the half-BPS sector in which the two-point function is diagonalised by the Schur 
polynomials Xr{ z ) [2]: 

(Xs{Z*)xr(Z)) oc 5 rs (1.3) 

The Schur polynomials were identified as gauge theory duals of giant gravitons [3], gen- 
eralizing the proposal of [4] which associated determinants and sub-determinants to giant 
gravitons expanding in the S 5 of the AdS§ x S 5 dual of J\f = 4 SYM. This is manifestation 
of the 'stringy exclusion principle' [5, 6]. 

x By 'matrix models in D spacetime dimensions' we include random matrix models (D = 0), matrix 
quantum mechanics (D = 1) or field theories. By 'Gaussian' we mean a quadratic Lagrangian. 
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The Schur polynomial basis allowed the identification of gauge theory duals for multi- 
ple giants as well as giants expanding in the AdS§ directions. Extensive evidence for the 
proposed map has been found, see for example the review [7]. It was also observed that 
Schur Polynomials were related to wavefunctions of free fermions which are the eigenval- 
ues of the matrix Z [2, 8]. The free fermions were subsequently identified as arising from 
supergravity solutions [9]. Non-renormalization theorems [10, 11] allow a comparison of 
zero coupling to strong coupling; the diagonalisation of the two-point function was a cru- 
cial step enabling this comparison. 

Motivated by the goal of understanding non-supersymmetric sectors of the AdS^/CFT^ 
gauge-string duality, [12] undertook the study of gauge invariant non-holomorphic func- 
tions of Z. Two-point functions of operators which are polynomials of degree m in Z 
and degree n in Z^ in the U(N) theory were diagonalised in terms of a Brauer basis 
which was constructed systematically using the representation theory of the Brauer al- 
gebra B]y(m,n), previously studied in [13, 14, 15, 16]. This non-holomorphic sector is 
arguably the simplest non-supersymmetric sector, yet contains many of the subtleties 
of multi-matrix combinatorics, since in general Z does not commute with its conjugate 
transpose Z^. 

The interest in a detailed study at zero Yang-Mills coupling is two-fold. Firstly, the 
strong form of the Maldacena conjecture [17] implies that there is a string theory dual to 
zero coupling Yang-Mills theory. At zero coupling, the sector containing only Z, Z^ forms 
a consistent truncation and we may ask whether there is a string dual of the quantum 
mechanics of the resulting complex matrix quantum mechanics. 

Secondly, while there is undoubtedly less control in comparing zero coupling with 
strong coupling and thus to semiclassical brane solutions and supergravity solutions it 
is possible that some qualitative features uncovered might, in appropriate large charge 
regimes, survive in the strong coupling limit. This line of reasoning is used, for example, 
in [18] where black hole entropy is investigated from the counting of gauge invariant 
operators at zero coupling. For some earlier works on complex matrix models, see for 
example [19, 20, 21, 22]. 

In this paper we always work at finite N. There is an important distinction between 

N > m + n and N < m + n . (1-4) 

The condition N > m + n may be read as a condition that N be larger than the lengths of 
operators one wishes to discuss; for example taking the planar limit N — > oo achieves this 
trivially. For fixed finite N, this is the regime in which lengths of operators are less than 
or equal to N. The opposite regime N < m + n is relevant to studies of heavy operators 
in Af = 4 super Yang-Mills such as conjectured duals of multi-branes and black holes in 
AdS 5 x S 5 . 
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The representation theory of Brauer algebras, and thus the construction of the Brauer 
basis, is well understood for N > m + n however there are interesting subtleties for 
N < m + n (see for example [23] ) . 

In single hermitian or unitary (or more generally, normal) matrix models, the unitary 
group action (1.2) is sufficient to diagonalise the matrix, leading to a free particle descrip- 
tion in terms of the N eigenvalues as we shall review. In an unrestricted complex matrix 
model, this is not sufficient to diagonalise Z\ a generic complex matrix may at best be 
put into triangular form. 

More precisely, using the Schur decomposition Z = UTU^ where T is upper triangular, 
the space gl(N; C) may be decomposed into a parameter space of inequivalent orbits M. n 
and the orbits of the U(N) action. A4n has real dimension N 2 + 1 and is a fibration over 
the symmetric product Sym N (C): 

M N 
I 

Sym N (C) = C N /S N (1.5) 

The eigenvalues are however coupled to the off-diagonal triangular entries and so cannot 
represent positions of free particles. 

We shall show that free particles arise in a non-trivial way by exploiting a map iden- 
tified in [12] between the k = sector (to be defined later) of the Brauer basis and a 
unitary matrix model, providing in turn a map to N free fermions on a circle. In this 
paper we will give evidence for the following conjecture: that these N free fermions of 
the k = sector can be constructed from degrees of freedom which are composed of 
eigenvalues as well as off-diagonal elements of the matrix Z. We also observe a different 
emergence of free particles in the m = n = k sector 2 . While we start with the gauge 
invariant sector of a Gaussian complex matrix model, which is a system of iV 2 particles 
constrained by the gauge invariance condition, the emergent particles are N free fermions 
without constraints. 

Much of our work has been carried out N = 2 as this allows explicit calculations. 
Many of our N = 2 results extend to general N as discussed in Section 5; in particular 
the key point of free particles emerging from a Matrix model from degrees of freedom 
beyond eigenvalues is valid for any N. 

1.2 Outline of paper 

The structure of the paper is as follows. Section 2 reviews the emergence of free particles 
from the eigenvalues of hermitian and unitary matrix models and explains the new features 

2 As this paper was being written up, we became aware of [24] which studies this sector and the 
associated free fermions using a matrix polar decomposition. 
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arising for a complex matrix model. Section 3 reviews the Schur decomposition, gauged 
matrix quantum mechanics and describes the orbits over M. n in general and over J\4. 2 m 
detail. In Section 3.5 we describe the ring of functions on M.2, which come from gauge 
invariant polynomials on gl(2,C). 

In Section 4 we describe the ring of Casimir operators studied in [25], we present 
computational results on the Brauer basis counting at N = 2 and state a conjecture for 
the complete solution. The conjectured counting can be elegantly described in terms of 
five integer labels and is the first main result of this paper. 

We derive explicit expressions for the Casimir operators as differential operators on 
M.2 and express the integer labels as functions of the Casimirs. We define free particle 
momentum operators and express these operators as functions of differential operators on 
M.2- Amongst these operators are the conjectured k = sector free fermion momenta on 
a circle. This is the second main result of this paper. 

Section 5 presents a conjecture that the k = sector is the kernel of a scale-invariant 
laplacian on Mn- We give expressions for a class of three-point functions of operators in 
the k = sector in terms of unitary matrix integrals, which provides further evidence for 
the conjectured equivalence to N free fermions on a circle. We extend some remarks on 
the counting of states at N = 2 to higher N. 

We also find free particle structures in the m = n = k sector, by observing that this 
sector consists of multi-traces of the combination Z^Z. We show that this sector may be 
identified with the kernel of a differential operator on Mn- 

Most of the material in Sections 3-5 refers to the properties related to invariants under 
the U (N) action (1.2) and as such is relevant to general complex matrix models with U (N) 
symmetry. Section 6 deals specifically with the matrix quantum mechanics obtained by 
dimensional reduction of SYM and draws a connection with Ginibre's D = matrix model. 
Using the higher conserved charges to define new Hamiltonians, and the expressions from 
Section 4, we find non-holomorphic generalizations of the Calogero-Sutherland models at 
special couplings. Technical details are presented in the Appendices. 



2 Review of free particles in Matrix models 

We briefly review examples of hermitian and unitary matrix models which arise in the 
context of string theory, in particular string theory in two dimensions. This review is not 
intended to be complete in any sense but rather to provide the reader with context for 
the current work. 



6 



2.1 Hermitian matrix quantum mechanics 

Let us consider the Gaussian hermitian matrix model defined by the Lagrangian 

£ = trQ$ 2 -i$ 2 ) (2.1) 

which is invariant under the U(N) action 

$ g$gi , g G U(N) . (2.2) 

We follow the treatment in [26, 27, 28] restricting attention to the theory with quadratic 
potential. The Hamiltonian of this model is 

H = tr(---^— + -<S> 2 ) • (2-3) 

Introducing the annihilation and creation operators 

A=±U + JL) A i = ±(*-JL) (2.4) 



(2.5) 



and using the usual convention for matrix indices 

d v 9 



we have [Aj, A' *] = 5 fe (^ and the Hamiltonian can be rewritten as 



iV 2 

i/ = tr(AU) + — . (2.6) 



The ground state has energy ^ and its wavefunction is 



2 

r tr<E> 



$ = ($|0) = e -a tr *' . (2.7) 

U(N) singlet excited states are obtained by acting on $o with U(N) invariant functions of 
A\ or by absorbing factors of a/2, multiplying by U(N) invariant functions of $. A basis 
for such functions is given by the Schur polynomials, which are polynomials of degree n 
labelled by a representation R of S n , 

XR (<S>)=J2xR(v)®%---®\ n , (2-8) 

cr£Sn 

where Xi?( cr ) is the character of a in the representation R. The associated wavefunction 

* R = XK($)e-^* 2 (2.9) 
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has energy ^- + n. 

A hermitian matrix $ may be decomposed as 

$ = UAlP, A = diag(X 1 ,...,X N ) , UeU(N) (2.10) 

under which we get 

tr($ 2 ) = tr(A 2 ) + tr[A, U^U] 2 . (2.11) 
The anti-hermitian matrix UW may be expanded in generators of U(N) as 



j<fc 

where Hi are the diagonal generators of the Cartan subalgebra, Tjk is the matrix M such 
that Mjk = Mkj = 1 and all other entries are 0, and is the matrix M such that 
Mij = —Mji = —i and all other entries are 0. This gives 

tv[A,U^U} 2 = J2(X i -X j ) 2 (al + ^ j ) 

i<j 

and so the Lagrangian becomes 

L = E fe 2 + + ^ 5> - A ;) 2 (4 + • ( 2 - 12 ) 

Under the transformation (2.2) the measure becomes 

£>$ = £>fi JJ dXiA 2 (X) (2.13) 

i 

where A(A) is the Vandermonde determinant rij<j(^« — the kinetic term for the 
eigenvalues becomes 

-l^A)i A2{ < t = -^a)?^ A(a) (2 - 14 » 

and so the Hamiltonian is 

» ~ 5E(-^(A) + A?)-ig^(^ + 4) .(, 15 ) 

Wavefunctions which are singlet under (2.2) are symmetric functions of the eigenvalues, 
Xsym{X). On these wavefunctions the Hamiltonian simplifies to 
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One may simplify further the analysis by defining the antisymmetric wavefunction 

tt'(A) = A(A) Xsym (A) (2.17) 

and the modified Hamiltonian 

fl, = A W fl zb = 5?(-^ + Aj ) (2 ' 18) 

which is a sum of one particle harmonic oscillator Hamiltonians. Then has eigenstates 
\1> / (A) with the same eigenvalues as H: 

HV(\) = EV(\) (2.19) 
=>H f V f (\) = Ety f (\) . (2.20) 

The ground state wavefunction of H* is 

^ = Ae^ tr *\ (2.21) 

excited states are given by Slater determinants 

H = detAjV^ tr * 2 = A(\)V R (U) (2.22) 

and so the U(N) singlet sector is equivalent to N non-interacting fermions in a harmonic 
oscillator potential, where the fermion energies are related to the integer row lengths 
Ti of R by 

Ei = n + {N-i) . (2.23) 
2.2 Unitary matrix quantum mechanics 

We next review the unitary matrix quantum mechanics which arises in the study of two- 
dimensional Yang-Mills, which is given by the Hamiltonian [29, 30]: 

H = tr(^) = J2 EUEa ( 2 - 24 ) 

where E a generates left rotations of U: 

E a = trt a U^- (2.25) 
oU 

The form of H means that acting on a wavefunction which is a matrix element of an 
irreducible representation R, 

WrMU) = D*(U) (2.26) 
9 



it measures the quadratic Casimir of the representation R, 

HMU) = C 2 (R)D*(U) . (2.27) 

Representations are classified by their characters, the Schur polynomials 

XR {U)=tTD R (U) (2.28) 

which form an orthonormal basis for wavefunctions invariant under the U(N) action 

U gUg\ g e U(N) . (2.29) 

This may be used to express any unitary matrix U as 

U = gDg\ D = diag{e ie \...,e u>N ), geU(N). (2.30) 

On functions invariant under (2.29), performing the change of variables (2.30) the Hamil- 
tonian becomes [29]: 

1 



-—A 
Ad6 2 



i L 

where denoting the eigenvalues by U{ = e ld 

-Ji = 



^N(N 2 - 1) (2.31) 



and where 

A(u) = Y[( Ui - Uj ) . (2.33) 

Absorbing A into the wavefunctions and the Hamiltonian, 

^ = , H t = AHl = 7J2 " ^ N ( N * ~ 1) ( 2 - 34 ) 

i 1 

the wavefunctions become antisymmetric under exchange of any pair 9^ -h- 9j. The one- 
particle wavefunctions with quantized momentum p are xf) p = e ip0 and the Slater determi- 
nants 

ipx = detu^ (2.35) 

i,3 

are eigenfunctions of Hf with energy E = Y2iPl ~ N(N 2 — 1)/12, so the sector of this 
theory invariant under (2.29) is equivalent to a theory of N free fermions on a circle. 
The ground state has fermion with momenta distributed symmetrically about n — 0, and 
energy zero, so the Fermi energy is np = — ^ and there are Fermi surfaces at ±np. The 
Slater determinants are related to the Schur polynomials via 

^ = A(u) XR (U) (2.36) 

where the momenta pi are related to the integer row lengths of R by 

Pi = U + (n F + 1 - i) . (2.37) 
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2.3 Complex matrix models 



Previous studies of complex matrix models have centred on models in which there is 
enough symmetry to diagonalise the matrix. This can be achieved by studying a normal 
matrix ([Z,Z<] = 0) with U(N) symmetry (see e.g. [31, 32]) 

Z -+ gZg^ , g e U(N) (2.38) 

or by studying an unrestricted complex Z with U(N) x U(N) symmetry (see e.g. [33]) 

Z^gZh ] , g,heU(N). (2.39) 

In this paper, motivated by gauge-gravity duality we study an unrestricted complex ma- 
trix Z with a single U(N) symmetry (2.38). This requires us to go beyond an eigenvlue 
description and take into account off-diagonal degrees of freedom. 

Due to the off-diagonal degrees of freedom we do not expect a straightforward transfor- 
mation to a description in terms of free particles for complex matrix models with unitary 
symmetry. Nevertheless, our investigations of the k — sector of the Brauer basis indicate 
that the free particles on a circle of the Unitary matrix model can be constructed from 
the degrees of freedom of the complex matrix model. These free particles on a circle are 
emergent degrees of freedom arising from eigenalues and off-diagonal elements constrained 
by equations which define the k = sector. 

In passing we note that the Gaussian complex matrix model with U(N) symmetry 
(2.38) may be written as a two-Hermitian matrix model [34] using 

X=±(Z + Z), Y = --{Z-Z) (2.40) 

where Z denotes complex conjugate of Z. Studies of the same model in terms of two 
hermitian matrices are done in [35, 36]. 

3 Orbits and parameter spaces 

The relation between gl(N,C), the space of complex matrices Z and the space M.n, of 
orbits under the adjoint action (1.2), is given by the Schur decomposition. 

3.1 Orbits and the structure of M. N 

Schur's decomposition (see e.g. [37]) allows one to write any complex matrix Z as 

Z = UTU ] (3.1) 
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where U G U (N) and T is upper triangular. It has been used previously in the context 
of the complex matrix model in [38, 39]. The eigenvalues Zi of Z become the diagonal 
entries (and hence the eigenvalues) of T. There are also off-diagonal elements ty for i < j. 
The equation (3.1) can be viewed as describing a map from the pair (U,T) to complex 
matrices. The map is onto, but not one-to-one. Pairs (U, T) and (e l9 U, T) describe the 
same Z. There is a U(1) N action 

U ->U'= UH, H = diag(e ie \...,e i9N ) 

T -+T'= H ] TH (3.2) 

which leaves Z unchanged. The diagonal e ld acts trivially on T but the ^(l)^ -1 part 
defined by 9j = mixes non-trivially with the angles in T. 

We can parameterize the coset U(N)/U(1) N using the variable L and decomposing 
U = LH (as for example in [40]) leading to 

Z = L(HTW)L ] = LTL ] (3.3) 

where T = HTW . It is also convenient to use the U(l) N ~ l part of (3.2) to set the N — 1 
entries on the superdiagonal of T (namely t/j+i) to be real, and to use (U,T). 

There is also the freedom, for fixed Z, to rearrange the eigenvalues in any order on the 
diagonal of T by altering U. This freedom exists because there is a Schur decomposition 
for each possible ordering of eigenvalues on the diagonal of T. Given 

Z = U^Ul = U 2 T 2 U\ (3.4) 

where T\ and T 2 have different orderings of diagonal entries, we have 

T 2 = (u\U^t(u\U^ = U 12 T X U\ 2 (3.5) 

where U\ 2 = U\U\. 

We have thus derived the construction mentioned in the introduction of M. n as a 
fibration over the symmetric product Sym N (C): 

M N 
I 

Sym N (C) = C N /S N . (3.6) 

The set of eigenvalues z 2 , . . . , z N of Z modulo permutations in forms the space 
Sym N (C). Local coordinates on the fibre of over Sym N (C) are obtained from the 
upper triangular elements ty, with i < j, appearing in T. 

Functions of degree n on M, N /S^ and natural inner products on the space of func- 
tions, which are expressible in terms of integrals, are organised by the symmetric group 
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S n . Since n can be arbitrarily large, we may say that Soo, defined as an inductive limit 
from finite symmetric groups (see e.g. [41]), is the symmetry organising the space of func- 
tions on M. n /Sn- In the case of A4n there is an infinite-dimensional underlying Brauer 
algebra constructed as a limit of finite algebras B^{m,n). 

The space of NxN complex matrices gl(N, C) consists of orbits generated by the U(N) 
action Z — > UZW. Due to the trivial U(l) action the real dimension of the parameter 
space of orbits M N is iV 2 + 1 = 2N 2 - (N 2 - 1). 

This suggests that the number of generators of ring of functions on M.^ should be 
N 2 + 1. This works in a straightforward way at N = 2, but in a nontrivial way at iV = 3. 
We will come back to this in Section 5. 

Local coordinates on A4.^ are given by and variables ty. At generic the orbits 
are topologically U (N) /U(l) = SU (N)/Z N . At ty = 0, the parameter space M. n becomes 
Sym N (C). The orbit is then generically SU(N)/U(1) N ~ 1 . Note that, when U(N) acts on 
its Lie algebra, the adjoint orbits are always Kahler (and hence even dimensional) [42]. 
This is no longer the case for orbits in the complexified Lie algebra gl(N, C). 



3.2 Differential Gauss's law 

Dimensional reduction of M = 4 SYM onto WL t x 5* 3 yields a U(N) gauged matrix quantum 
mechanics involving a complex matrix Z(t) in the adjoint coupled to a gauge field Ao(t). 
The action takes the form 



S = J dt tr ^D Z(D Zy - ZZ^j 



(3.7) 



where D Z = d Z + i[A , Z]. 

Using the above change of variables (3.1), (3.3) we may derive an expression for the 
1-form on gl(N; C) 

dZ = U(dT + [u, T])W = L(df + [V, f ])L f (3.8) 

where u = U^dU and V = L^dL. This allows us to write the line element tv(dZdZ') in 
terms of the structure constants of the Lie algebra, with a choice of decomposition into 
coset and sub-algebra. 

As an aside, it is interesting to note that, as a consequence of (3.8) we can write a 
quantum mechanics theory with U (N) global symmetry 

S = J dt tr (d Q Z 8 Z^ = J dt tr (dot + [V, f ]) (d f f + [V, ft]) (3.9) 

as a quantum mechanics with gauged U(1) N symmetry and charged matter fields T where 
the one form V on the coset couples as a gauge field. The gauge symmetry is T — > hTh~ x 
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and V(y) —> hVhr x + hdoh' 1 , under which (<9oT + [V, T]) transforms covariantly and the 
action is invariant. 

We next review remarks contained in [25] and introduce notation we shall use later. 
A convenient gauge fixing choice is to set A = 0. The equation of motion for A must 
still be imposed, leading to Gauss's Law: 

Z*Z + ZZ^ - ZZ* - Z X Z = . (3.10) 

Upon canonical quantization this leads to the differential form of Gauss's Law, which can 
be written as 

G = Gi + G 2 + G 3 + G 4 = (3.11) 

where Gi are defined as: 

< G '>'< = z ^ W)] W> - *> G§)" 

^ = {mj t < G '« = " z * W„ (3 ' 12) 

and we use the usual convention for matrix indices given in (2.5). Note that in G\ and 
G2 the ordering of indices is that of usual matrix multiplication, while for G3 and G4 the 
opposite is the case. The Gi correspond respectively to each of the terms in (3.10). The 
operator G is the infinitesimal generator of the adjoint action 

Z -»• UZU\ Z ] -> Utfrf (3.13) 

and invariance under this action restricts gauge invariant operators to be products of 
traces of the matrices Z and Z\ 



3.3 Geometry of M.^'- coordinates 

In this section and in Section 4 we perform explicit calculations at N = 2. The motivation 
for considering small values of N is to perform explicit calculations which shed light on 
the harder (and more interesting) task of obtaining results at arbitrary finite N, a task 
we return to in Section 5. 

We start from the Schur decomposition as discussed in Section 3.1, 

Z = UTU^ = LTl) . (3.14) 

In the N = 2 case U{2)/U{1) ^ SU(2)/Z 2 = SO (3). We can specify explicit coordinates 

n ( cosfe^) sinfe^)\ 

^_ sin | e -f(^) cos|e-3(^y 1 } 

T = (* . (3 , 6 ) 



14 



The angles 0,cj),ip are the Euler angles of SU(2)/Z 2 = 5*0(3). With these coordinates L 
and T take the form 



cos^ea^ sin ^ e 2 s 
-sinle - ^ cos|e~; 



z 1 t e # 
z 2 



(3.17) 
(3.18) 



The ranges of the coordinates are 



Zx, z 2 eC, < t < oo, (3.19) 
< 6 < 7T, < < 2vr, < ip < 2vr. (3.20) 

The Jacobian for the change of variables from Zij to those above is 

J —\z 1 — z 2 \ 2 t sin6> (3.21) 

and so we have 

J W^^^^v = J dzidz 2 dtotodU \zi — z 2 \ 2 . (3.22) 

Note the factor of to here which is analogous to the J rdr one gets when using plane polar 
coordinates. Here dll is the Haar measure on SU(2) which we integrate out and normalise 
to 1 in the definition of the measure. The implication of the measure is that the region 
to = 0, where the orbit structure changes compared to that at to 7^ 0, has measure zero. 
Likewise the collision of points z\ = z 2 in Sym N (C) has measure zero. 
The invariant line element on gl(2, C) is given by 

ds 2 = tidZdZ ] . (3.23) 

We introduce the notation 

u = U-'dU = ( Un Ul2 ) , (3.24) 

and using = — u we expand dZ = U (dT + [u>,T])W. 
The line element is then expressible as 

ds 2 = tr (dT+ [u,T])(dT^ + [u,T*\) (3.25) 
= \dzi + t u 12 \ 2 + \dz 2 - t ui 2 \ 2 

+ \dt + 2t wii - (zi - z 2 )uu\ 2 + \(zi - z 2 )u 12 \ 2 . (3.26) 

Using the Cartan one-forms Ui on SU(2) (see e.g. [43]), 

u = U~ l dU = -UiT i7 Tj = l -<Jj, (3.27) 

15 



one may read off the metric on the orbit; we shall do this in the next section. 

As an aside, we note that {7 12 defined below 3.5 is not a standard permutation matrix 
in U (N) (the reader may check that the standard permutation matrices in U(N) do not 
preserve the triangular form). For concreteness we now exhibit this at iV = 2. Consider 
the two matrices 

71 = (* <°) (3.28) 
T 2 = (* j (3.29) 

where we have chosen t e KL 

Defining D = y/t^ + \z~\ — z 2 \ 2 , we then have T 2 = U\iT\U\ 2 with 

= ~V 2 ')' (3 ' 30) 

Clearly this is not the standard permutation matrix (JJ), but it performs the permu- 
tation transformation £i -H- z 2 while preserving the triangular structure. For N > 2 
the analogous transformation does not just permute the Zi entries but transforms the tij 
nontrivially. 



3.4 Differential Gauss's law and orbits at N = 2 

Using a change of variables, one may express the Gauss Law operator G (3.11-3.12) in 
the coordinates defined in (3.15-3.16). This results in the following form of the Gauss's 
Law operator: 

r _( ie^C-^-tcot^ + icsc^) \ 

G -\ie-^ e -icotel + icsce±) i± )■ ^ 

This must vanish on gauge invariant wavef unctions, which must therefore be functions 
only of Zi,t . We will show in Section 3.5 that the ring of gauge invariant polynomials 
has five generators. 

The Gauss's Law reduces the 8D space gl(2,C) to the 5D space parametrized by 
(zi,z 2 ,to). We shall find it convenient to define 

z c = zi + z 2 , z = zi - z 2 . (3.32) 

As we have seen, we can exchange Zi,z 2 while leaving to invariant; this means mapping 
z — » —z, and so the space of inequivalent orbits is 

M 2 = Cx (C/Z 2 ) x E+. (3.33) 
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From the metric (3.26) expressed in terms of Zj, to we see that the nature of the orbits 
changes as we move in the space (C/Z 2 ) x R + . The centre of mass coordinate z c does not 
affect the nature of the orbits and so we restrict our attention to a Z 2 quotient of the z, to 
space. Let us define 

X = (C/Z 2 ) x R + 

= X UX 1 UX 2 UX 3 (3.34) 

where X is the region in (z, to) space where to > 0, Re(z) > 0, and the subregions Xi are 
defined as follows: 

• X is the subregion t > 0, z ^ 

• X\ is the subregion t > 0, z = 

• X 2 is the subregion t = 0, z ^ 

• X 3 is the point to = 0, z = 0. 

The metric on the gauge orbit is determined by fixing z^ to in (3.26). On X Q and 
Xi the orbit is topologically SO(3); the metric is complicated in general but on X\ it 
qualitatively resembles the round three-sphere metric. On X 2 the orbit is a round S 2 , 
while on X 3 the orbit is a point. This completes the global description of the parameter 
space and the orbits. Note that on Xo the metric is regular but on X±, X 2 and X$, the 
determinant of the metric is zero. 

3.5 The algebra of functions on M.2 

The algebra of functions on M. # is generated by single trace polynomials in Z, Z* . In the 
N — > oo limit any word in the two letters Z, Z', up to cyclic permutations, corresponds 
to a single-trace gauge-invariant function and hence to a function on .Moo- At finite N, 
traces of long words can be expressed in terms of products of traces of shorter words and 
so the ring of gauge invariant functions has a finite set of generators. 

In [12] this truncation of the generators was discussed in terms of degenerations of 
Brauer algebra projectors. Here we investigate these finite N truncations in detail at 
N = 2 and find that it suffices to apply the Cayley-Hamilton theorem to obtain the 
necessary relations. 

The Cayley-Hamilton theorem states that a matrix satisfies its own characteristic 
polynomial. At N = 2 this means that 

Z 2 - (trZ)Z+ (detZ)h = 0. (3.35) 

Taking the trace of this equation gives a relation between tr Z 2 , tr Z and det Z, only two 
of which are thus algebraically independent as polynomials in the matrix entries. We 
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choose tr Z 2 and trZ to be independent, and write 

detZ = - [trZtrZ - trZ 2 ] . (3.36) 
We also have the corresponding equation for Z\ 

We claim that the ring of multi-trace GIOs in Z, at N = 2 is the polynomial ring 
generated by the set 

£={trZ, trZ 2 , trZ f , trZ f2 , trZZ f }. (3.37) 

In order to prove this, it is enough to show that all other single trace operators are 
algebraically dependent on the operators above. 

We prove this in an inductive fashion. Let W to denote any matrix word made from 
Z and Z\ e.g. W = ZZZ^Z. Multiply (3.35) by W and take the trace. This yields the 
relation 

tr(Z 2 W) - (tr Z) tr(ZW) + - [tr Z tr Z - tr Z 2 ] tr W = 0. (3.38) 

This shows that tr(Z 2 W) is algebraically dependent on ti(ZW), tr W and the operators in 
B, and similarly, ti(Z^ 2 W) is algebraically dependent on tr(Z^W), trW and the operators 
in B. 

Replacing Z by ZZ* in (3.35) and using det ZZ* = det Z det Z^ gives 

tr(ZZ^) 2 = (ti ZZ ] ) 2 - - [trZtrZ-trZ 2 ] [tr tr Z f - tr Z t2 ] . (3.39) 

This shows us that tr(ZZ^) 2 is algebraically dependent on the operators in the set B. 
Similarly, for any word W 2 of length at least two, tr W% is algebraically dependent on 
tr W2 and the operators in the set B. 

We conclude that a single trace operator consisting of the trace of a word made from Z 
and Zt is algebraically dependent on single trace operators of shorter length iff it contains 
one of the following combinations as part of the word: 

Z 2 W, Z ]2 W, or Wl (3.40) 

where as above W stands for any (non-zero length) word in Z and Z\ and W2 stands for 
such a word of length at least two. 

Iterating the above results, a single trace operator containing one of the combinations 
in (3.40) can be expressed as sums of products of shorter and shorter single trace operators 
until it is expressed as a sum of products of single trace operators containing none of the 
combinations in (3.40). A maximal set of algebraically independent operators is therefore 
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given by those single trace operators which do not contain any of the expressions in (3.40). 
As claimed this is the set B. 

It is worth remarking that we start with a description of the space gl(2, C) in terms 
of polynomials in zi, z 2 , t , 6, <j), ip. The differential Gauss Law (3.31) removes the angular 
variables leaving the ring of polynomials in the remaining variables, which we denote 

(z 1 ,z 2 ,z 1 ,z 2 ,t ). (3.41) 

Invariance under large gauge transformations reduces the ring of gauge invariant polyno- 
mials to the polynomial ring generated by B. Recalling the definitions z c = Z\ + z 2 , z = 
Z\ — z 2 and defining 

Z = z\ Z = z\ T = tg + y, (3.42) 
the ring of gauge invariant polynomials is equivalently the polynomial ring 

(z c ,z c ,Z,Z,T ). (3.43) 

This is analogous to U(N) gauged Hermitian matrix quantum mechanics where the dif- 
ferential Gauss Law reduces to polynomials in the eigenvalues 

(x u x 2 , . . . ,x N ) (3.44) 

and invariance under the Sn residual Weyl transformations reduces the gauge invariant 
polynomials to symmetric polynomials in xi,x 2 , ■ ■ ■ , xjv, equivalently polynomials in the 
variables 

((a* + x 2 + ■ ■ ■ + x N ), {x\ + x\ + ■ ■ ■ + x%), . . . , (xf + x% + ■ ■ ■ + a#)). (3.45) 

In the hermitian case, we are going from a ring to a sub-ring, which corresponds to going 
from the space M N to its quotient space WL n /Sn- In our model, we are going from the 
ring (3.41) to the sub-ring (3.43), and correspondingly from the M 4 x M + = C 2 x IR + 
parametrized by the five coordinates Zi,t to M. 2 . Because of the off-diagonal degrees of 
freedom, Ai 2 is not a straightforward quotient of M 4 x M + . 

A full investigation of finite N relations for N > 2 is left for the future. We expect 
it will be useful to combine the Cayley-Hamilton approach with the the vanishing of the 
Brauer projectors, such as in equation (8.16) of [12]. 

4 Free particle structures and counting on M.2 

The remainder of the paper involves the Brauer basis for complex matrix models con- 
structed in [12]. A Brauer basis operator is a linear combination of multi-trace operators; 
it is a polynomial of degree m in Z and degree n in A brief review of the essen- 
tial properties of the basis and some simple examples are given in Appendix A, where 
references to existing literature are also given. 
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A Brauer basis operator is written as 



(4.1) 



where a and are Young diagrams with m and n boxes and 7 = (/c,7 + ,7_) where k is 
an integer in the range < k < min(m,n) and 7+, 7- are Young diagrams with m — k 
and n — k boxes respectively. For a more complete explanation of the labels please see 
Appendix A. 

As an example for the reader to bear in mind, when (m,n) = (1,1), suppressing 
non-essential labels the Brauer basis is 



Since we discuss in particular the label k throughout the rest of this paper, the following 
comment will be useful to the uninitiated reader. In the construction of the Brauer basis, a 
term with a single 'ZZ' J inside the same trace, such as tiZZZ\ involves a single 'Brauer 
contraction'. Terms such as tr ZZ^ tr ZZ^Z^ or tr ZZ^ZZZ^ involve two such Brauer 
contractions, etc. 

The label k is related to the number of contractions as follows. If one writes a Brauer 
basis operator as a sum of terms in order of increasing contractions, as the two operators 
above are written, an operator with label k begins with a term involving k Brauer con- 
tractions. We have not proved this, but we believe it to be true from all the examples we 
know. Thus the leading term in a k = operator is the product of a purely holomorphic 
operator and a purely anti-holomorphic operator, while all terms in a m = n = k operator 
involve k contractions. 

The first result of this section is a conjecture for the solution to the N = 2 counting 
of the operators of the Brauer basis, for which which we provide numerical evidence. 

The second result will be to find evidence of a 'free fermions on a circle' structure in 
the k = sector. This generalizes to any N, as discussed in Section 5.1. We show that 
the Brauer basis at A^ = 2 can be neatly expressed in terms of five integers and observe 
the correspondence between states in the k = sector and two free fermions on a circle. 
In these developments a crucial role is played by the structure of the ring of Casimirs. 

The third result in this section is to show that the momenta of the free fermions of the 
k = sector can be constructed from differential operators in variables which include both 
eigenvalues and off-diagonal elements of Z. This leads us to observe that the complex 
matrix model contains free fermions arising in a novel way, different from the way they 
arise in hermitian or unitary models. 



trZtr Z ] - ^-tiZZ ] 



(4.2) 




(4.3) 
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4.1 Casimir operators and a ring of degree-preserving differen- 
tial operators 

The differential operators introduced in equation (3.12) were studied in [25] as generalized 
Casimirs commuting with the scaling operator for Z, Z', which is the Hamiltonian for zero 
coupling SYM. This ring is analogous to the ring generated by B in Section 3.5; at N = 2 
the generating set is 

X?={trG 2 , trGl trG 3 , txG\, trG 2 G 3 } (4.4) 
where G 2 ,G 3 were defined in (3.12) 

m = Z'J^j ( G3 ).=-W^J (4.5) 



Defining 



we introduce the Hamiltonians 



G L = G 2 + G 3 , (4.6) 



iZi=trG 2 H 2 = XiG 2 2 



H x = tr G 3 H 2 = tr G£ H L —trG L . (4.7) 

Each of these operators commutes with the scaling operator for Z and Z\ which is 
H = Hi + Hi. The operators in T> generate a ring of commuting Hamiltonians related 
to the integrability of the system. We have defined Hl for later convenience; its name 
derives from the fact that the operator G 2 + G 3 is the infinitesimal generator of the left 
action of U(N) [25]: 

Z -»■ UZ, Z ] -»■ Z ] U ] . (4.8) 
It was shown in [25] that the five operators defined in (4.7), 

H a = {h 1 , Hi, H 2 , H 2 , H l } (4.9) 

measure respectively the Casimirs 

'Ci(a), d(/3), C 2 (a), C 2 (/3), C 2 (j)\ . (4.10) 



Generalized Casimir operators such as tr(G2G 3 ) were investigated in [25] and were shown 
to be sensitive to the labels i,j in (4.1). Since the matrix elements of G 2 and G 3 commute, 
we may regard G 2 and G 3 as matrices of c-numbers and apply the Cayley-Hamilton 
theorem as in Section 3.5 to show that the set T> is a maximal algebraically independent 
set of degree-preserving gauge invariant differential operators. 
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This observation implies that the generalized Casimir operators such as tvG^Gs do 
not yield independent information about the wavefunctions at N = 2, i.e. that all the 
information in the labels {a, /3,7, i,j} is in fact contained only in {a,/3, 7}. We can 
interpret this fact in terms of Brauer algebra representation theory as follows. 

In the restriction of an irreducible representation 7 of the Brauer algebra to the rep- 
resentation A = (a, (3) of C[S m x S n ], there enters an integer multiplicity M^ N defined 
by 

yB N (m,n) _ M]^ N Y^ SmXSn ^ (4 ]7) 

A 

For large N, i.e. m + n < N, we denote this multiplicity by M\ or M^p and using 5 h k 
to denote that 5 is a partition of k we have the formula [14] 

M l = = EE 9(l+A*)9<n-AP) (4-12) 

<5hfc 5 

where #(7+, 5; a) is a Littlewood-Richardson coefficient. 

As reviewed in Appendix A the indices i, j on a Brauer operator range over the values 
{1, . . . , M^ N }, and so the redundancy of the i, j labels at iV = 2 means that M"^ N ~ 2 is 
either or 1 for all 7, A. A direct proof of this by using the finite N constraints on the 
states of the Brauer representation in [15] would be interesting to obtain. At this point 
we will take a more pragmatic perspective, assume it is true, and will find that it leads 
to a consistent counting of states of the complex matrix model at N = 2. 



4.2 Counting of states at iV = 2 and Brauer basis labels 

The ring of gauge invariant operators at N = 2 is generated by five single trace operators 
(3.37). Hence the number of linearly independent multi-trace operators Q^f 2 (m,n) for 
fixed (m, n) is counted by the generating function 



\ r 

(1 - x){l - y)(l - x 2 )(l - y 2 )(l - xy) ^ 



Q N m f{ mi n)x m y n . (4.13) 



rnji 



This is the Plethystic Exponential [44, 45] of the single trace generating function 

J2QT 2 ( m > n )% m y n = l+x + y + x 2 + y 2 + xy (4.14) 

m,n 

derived from the independent single traces in the basis B (3.37). 

Having found the N = 2 counting of multi-traces, we can express it in terms of 
constraints on the large N Brauer counting. The obvious constraint £7(7+) + Ci(7_) < 2 
is not sufficient. We have argued above that the multiplicities M^~ 2 are either or 1. 
We first set 



otherwise 
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where Mj^ is given by (4.12). Having done this we also find it necessary to impose extra 
constraints on the labels a,/3 for agreement with (4.13). 

The constraints on a,/3 are as follows. Denoting the length of the p th column of a 
Young diagram R by c p (R), we constrain: 

1. ci(a) + Ci(/3) < N + k 

2. [ci(a) + Cl (/3)] + [c 2 (a) + c 2 (/3)] < 2iV + k 

and in general for each p = 1,2, . . . , min(m, n) , constrain 

v 

Y^{c r {a) + c r {l3))<pN + k. (4.16) 

r=l 

We have used SAGE and Mathematica to enumerate all possible Brauer basis operators 
subject to the constraint (4.16) and to compare with the Trace basis generating function. 
The two agree up to (m,n) = (15, 15) which is the practical limit for a desktop computer. 
This conjecture generalizes the 'Non-chiral Stringy Exclusion Principle' introduced in 
[12]. This counting of operators at iV = 2 implies a result for the reduction multiplicities 
M^ N ~ 2 , namely that 

M r,N=2 f 1 if > and (4.16) holds ^ 

a,/3 I otherwise 

We will re-state this result after simplifying the condition (4.16). 

We can use the fact that the Brauer basis diagonalises the five Casimirs %a to explicitly 
enumerate the independent Brauer operators at N = 2. Fixing (m, n) we pick a basis 
of multi-traces. Acting with the explicit form of the Casimirs (4.41), we can construct 
linear combinations which are eigenstates of the Casimirs. Because the eigenvalues of the 
five Casimirs determine the labels a, /3, 7 uniquely, we can also read off the labels of the 
allowed operators. We have carried out this procedure for selected values of (m, n) up to 
(m,n) = (4,3). 



4.3 The Brauer basis labels at N = 2 in terms of five integers 

In Section 3.5 we described the states of the N = 2 theory as generated by a finite set 
of traces. In this section we will obtain the description in terms of the Brauer basis for 
multi-traces. For general N, we give a review of the Brauer basis states in Section A. For 
ease of notation we denote = rj(a) and fj = Vi(j3). 
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We can choose different sets of five integers to parameterise the states, such as 

r 1 ,r 2 ,f 1 ,f 2 ,rj (4.18) 
rj, r 2 , k, r±, f\ (4.19) 
r1,r 2 ,k,r\,r 2 . (4.20) 

We will show that each of the above sets of five integers determines a state uniquely, and 
we will give the constraints on the integers. 

A state is determined uniquely at N = 2 by a, (3,j, containing the set of integers 

{n,r 2 ; fi,f 2 ; k,rj,r^}. (4.21) 

From the Brauer algebra representation theory briefly reviewed in Appendix A, we have 
the following relations : 

= m, = w, (4.22) 

i i 

^ n{j + ) =m-k, r i{l-) = n — k. (4.23) 

i i 

Using the relationship between 7^(7), ^(7+) and ^(7-) we have 

53^(7) = ^2 r i(l+) -^2 r i(l~) = rn-n (4.24) 

i i i 

which at N = 2 reads 

r i + r 2 = m ~ n - (4.25) 
Adding the two expressions in (4.23) we find that 

^Vi(l)\ = Y. ri ^ + Ys ri ^^ = m + n ~ 2k ( 4 - 26 ) 

i i i 

which at N = 2 gives 

k = -( m + n-\rj\-\r^\). (4.27) 
We now show that each of (4.18)-(4.20) are enough to determine the state via (4.21): 

1. Starting from the five integers in (4.18), we deduce m,n from (4.22), from (4.25) 
and k from (4.27). 

2. Starting from (4.19) we read off ^(7+) and rj(7_) by inspecting whether and rj 
are positive or negative. We then deduce m and n from (4.23) and r 2 and f 2 from 
(4.22). 

3. Starting from (4.20) we proceed as in point 2 above. 

This shows that each of the three sets of five integers identified are sufficient to identify 
any state. 
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4.3.1 N = 2 constraints in terms of five integers 

Let us consider the case where k is one of our five integers. We rewrite the N = 2 
constraint (4.16) as a lower bound on k: 

v 

k > ^J(c,.(a) + c r (f3)) — 2p for each p — 1, . . . min(m, n). (4.28) 

r=l 

Note that as p increases the lower bound on k gets stronger only when 

cp(a) + c p ((3) > 2. (4.29) 

Before presenting a general expression for the lower bound on k we examine in detail the 

case 

< r 2 < f 2 < ri < fi. (4.30) 

We observe that 

• For 1 < p < r 2 we have c p (a) + c p (/3) = 4 

• For r 2 < p < ?2 we have c p (a) + c p (/3) = 3 

• For p > f 2 we have c p (a) + c p (/3) < 2 

The strongest lower bound on k is therefore at p = f 2 where we have 

k > 4r 2 + 3(f 2 - r 2 ) - 2f 2 
=^ A; > r 2 + f 2 . (4.31) 

Proceeding similarly we find a general expression for the lower bound on k. For simplicity, 
wlog suppose r 2 < f 2 . There are three cases to consider: 

1. r 2 < ri < f 2 < f\ k > ri + r 2 

2. r 2 < f 2 < r! < f\ k > r 2 + f 2 

3. r 2 < f 2 < r 7 ! < r! ^> /c > r 2 + f 2 . 
Combining these we obtain the lower bound 



k > min(r 2 ,r 2 ) + min ( min(ri,ri), max(r 2 ,r 2 ) ) (4.32 



which is equivalent to (4.16). We can also express the constraint (4.32) in terms of the 
five integers in (4.18) by substituting for k from (4.27) to find 

— (m + n — \rj\ — \m — n — rj\) > min(r 2 , f 2 ) + min ( min(ri, fi), max(r 2 , f 2 ) ) . (4.33) 
We can now re-state the result (4.17) for the N = 2 reduction multiplicities: 

(4.34) 



M r, N =2 ] x if K,p > and ( 4 - 32 ) holds 



[ q,/3 



otherwise. 
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4.4 The Casimirs as differential operators in to 

In this section we express the Casimir operators / Hamiltonians from Section 4.1 as 
differential operators on Ai 2 . 

Below are calculated expressions in the coordinates Zj,£o f° r the Hamiltonians defined 
in (4.7). For convenience define 

d d 
L x = zi— Li = z 1 — (4.35) 

OZi OZi 

u = z *ik u= ~ Z2 h Lt = t iw (436) 

and recall the notation z c = z\ + z 2 , z = z\ — z 2 . 

Recalling the definition Gl = G2 + G3 from above (4.7), we find the following expres- 
sions: 

H 1 = tiG 2 = Li + L 2 + L t (4.37) 
H x = -trG 3 = U + Lz + Lt (4.38) 

H 2 = trGj = Ll + Ll+^-^f)^ 

+ -{z 1 L 1 -z 2 L 2 )L t + ^{L 1 -L 2 )+L t (4.39) 



H 3 = tiG 2 3 = tr(G|) (4.40) 

H L = trG\ = (L x - LO 2 + (L 2 - L 2 f + ^ (L x - L 2 ) + ^ (Z x - L 2 ) 

z z v ' 



2 f . 9 - ^ , 1 

r o 



t (Li - L 2 )(Li - L 2 ) + -2(^1^1 - 2^2) L t 

\Z I In 



-(zizi - z 2 z 2 ) [(Li - L 2 ) + [L x - L 2 )] L t - (z 1 z 1 + z 2 z 2 )L t j 



(4.41) 

Some useful formulae in doing these calculations are now given. Recall from (3.8) the 
definition V = L^dL and the expression 

dZ = L (df + [V,f]\ Ll (4.42) 

Defining 

dX = df+[V,f\ and (G 2 )) = T p (JVj (4.43) 
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one may derive 



d_ 

dZ 



dZ) = VpdX^L^ 

p 



(G 2 )) = L\p {G 2 )\. (4.44) 
The computation of (G 2 ) p shows that it contains angular derivatives. When we calculate 

trG* = UplPj (G 2 ) p q L i r L s j (G 2 ) r s (4.45) 

it is important not to neglect the terms obtained from the action of these angular deriva- 
tives from {G 2 ) p q on L l r L s -. 

4.4.1 The Casimirs as operators on polynomial rings 

We observed in equation (3.43) that the multi-trace operators built from Z, Z' form a 
polynomial ring whose generators we may take to be 

z c , z c , Z = z 2 , Z = z 2 , T = ^ + y. (4.46) 

The above differential operators H 2 , H 3 , map polynomials in these variables to poly- 
nomials. Changing variables to these generators makes this manifest: 

H 2 = 2l(l + ^+±L c (L c + 3) + L (Lo+1) + ^L + £(2L-1)L 

+ ^L C (L C - 1) + ^(z 2 - Z)L (L - 1) + 2 + ^) LL + 2L L C + ALL C 



(4.47) 

where L = Z-J^, Lq = Tq-^t and L c = z c -J=^-. H3 is obtained by complex conjugation and 
the same exercise can also be done for Hi to illustrate that they are operators that map 
polynomials to polynomials. 

4.5 Eigenvalues of the Casimir operators 

As reviewed in Section 2, a Young diagram R with non-negative row lengths labels 
energies £j of N fermions in a one- dimensional harmonic oscillator potential, given by 

£i = n + (N- i) (4.48) 

and a Young diagram (iV-staircase) R with arbitrary integer rj labels momenta pi of N 
free fermions on a circle given in terms of the Fermi energy rip = by 

Pi = n + (n F + 1 - i) . (4.49) 

In this section we review the fact that the values of 
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• the N independent U(N) Casimirs Ci(R) of the representation R 

• the N row lengths r i; and 

• the N corresponding fermion momenta Pi are equivalent data. 

The same remark holds for non-negative with replaced by £j. 

In Section 4.1 we introduced differential operators studied in [25] which when acting on 
a Brauer basis function O^AZ, Z^) measure the quadratic Casimir of the Young diagrams 
a,/3,7. Given a U(N) Young diagram R, its linear and quadratic Casimirs are 

dOR) = J2 r * = n ( 4 - 5 °) 

C 2 (R) = nN + ^V i (r i -2z + l). (4.51) 

i 

Using the definition of pi (4.49) we can write C 2 as 

N N 

C2{R) = E^-T2 (Ar2_1) (452) 

i=l 

which agrees with (2.34). Using the definition of Si (4.48) we can also write C 2 as 

N N 
C 2 (R) = ^^-(iV-l)n--(iV-l)(2iV-l). (4.53) 

i=l 

For general N, knowledge of the values of the independent Casimir invariants Ci de- 
termine the values of the power sum symmetric polynomials 

V a = p a 1 +P a 2 + +P% (4.54) 

which in turn for a — 1, . . . , N enables us to solve for pi or respectively Si (see e.g. [46]). 

We now demonstrate this in the N = 2 theory. The free fermions on a circle have 
ground state with energy pi — |, p 2 = — | and in general we have 

1 1 . 

Pi = n + -, p 2 = r 2 - - . (4.55) 

Setting N = 2 in (4.51) gives 

C 2 = n{ ri + 1) + r 2 (r 2 - 1) (4.56) 
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and so we may express C\ and C 2 in terms of pi as 

Ci = pi+ p 2 

C 2 = pl + pl-±. (4.57) 
The resulting quadratic equations for pi in terms of C\ and C 2 have solution 

C x , /c^ C? ' 1 

Pl = T + VT-T + 4 

p 2 = Ci-px. (4.58) 

4.6 The jfc = sector 

In the k = sector 7 = (0, a, /3) so operators are labelled simply by a and /3 which are 
representations of S m and S n respectively. To connect with the notation of the unitary 
matrix model, we write a = R and (3 = S. If S = 0, then the k — operator is the 
holomorphic Schur polynomial corresponding to the representation R: 

O k R : d °(Z,Z^ = XR (Z) . (4.59) 

If R = 0, then the = operator is the anti-holomorphic Schur polynomial corresponding 
to the representation S: 

0$f(Z,&)= Xs (tf) (4.60) 

and if both a and f3 are nontrivial, the leading order term in the expansion of O k=0 begins 
with the product of the holomorphic and antiholomorphic Schur polynomials: 

O k =§{Z, tf) = X r(Z)xs(&) + • • • , (4.61) 

where the dots denote terms with at least one ZZ^ inside a trace as discussed at the start 
of Section 4. 

There is an isomorphism between the k = sector and the states of the Unitary matrix 
model [12]: 

6${Z,Zt)<-^ X RS{V) (4-62) 
which is obtained by replacing Z with a unitary matrix: 

O k R f(U,U^) = d R d sXR s(U) . (4.63) 

The two point functions of both sets of operators are diagonal; up to a choice of normal- 
isation, 

{0 R k -:\Z,tf)\O k =!,{Z,rf)) = (xUu)\ X r>s>(U)) = 5 RR ,5 § s, (4.64) 
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and the reader familiar with the 'coupled characters' studied in two-dimensional Yang- 
Mills will notice that the structure of (4.61) is of the same form as the coupled character 
Xrs- The k = states are thus isomorphic to the states of N free fermions on a circle via 
the map given in the same section. 

At N = 2, the label 7c as defined in (A. 4) may have at most two rows, r\,r\ and 
so the integers (k = 0, rj, r 2 ) are enough to specify an operator. The list of all N = 2 
operators for given (m,n) in Appendix B shows that: 

• If rj > 0, r\ > 0, then /3 = and we have a holomorphic Schur polynomial. 

• If rj < 0, r\ < then a = and we have an antiholomorphic Schur polynomial. 

• If r\ > 0, r\ < then the operator is of the form (4.61). At N = 2 there is a unique 
such operator. 

Since row lengths and fermion momenta are equivalent data in specifying a state, the 
above constraints may be rewritten in terms of fermion momenta p]. In the next section, 
we will see how the momenta of these fermions can be expressed in terms of differential 
operators in Zi,t . 

4.7 Free particle momenta as functions of differential operators 

As noted in (4.9), when applied to an N = 2 Brauer basis operator O^p, the differential 
operators 

Ha={h 1 , H u H 2 , H 2 , H l } (4.65) 
measure the values of the Casimirs 

C A = {Ci(a), Cxi/3), C 2 (a), C 2 (/3), C 2 ( 7 )} (4.66) 

respectively. We also have the fact that 6*1(7) is measured by Hi — Hi. We define fermion 
momentum operators 

Pa = {pi, P2, Pi, Pa, Pi, P2} ( 4 -67) 

whose eigenvalues are pi, p 2 , pi, pi, p\, p 2 respectively. We now repeatedly apply (4.58) 
to each of a, (3, 7 in turn which enables us to derive expressions for these operators in terms 
of the basic gauge invariant operators Ha- 

Applying (4.58) to the label a and promoting to an operator equation we obtain 

Hi [W 2 Hi I 
p 2 = Hi-pi. (4.68) 
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Applying (4.58) to the label (3 we obtain analogous expressions for p±,p2 in terms of 

Applying (4.58) to the label 7, promoting to an operator equation and defining d = 
Hx — Hi we obtain 

d [h~ l J 2 1 
ft = 2 + V^~T + 4 

Pl = d-pj (4.69) 

As noted in Section 4.6, in the k = sector a state is specified simply by the values 
of the row lengths rj, rj, or equivalently by the values of the fermion momenta p\-,p\ and 
so we now identify pj,p2 as formal expressions for the momenta of the k = fermions on 
a circle. We shall extend this result to arbitrary N in the next section. 

Comparing to the explicit expressions for T-La obtained in Section 4.4, we see that 
these fermion momenta are functions of differential operators in both the eigenvalues 
and the off-diagonal element t - I n hermitian matrix models and unitary matrix models, 
the emergent fermions are the eigenvalues of the relevant matrix. Here, however, the 
k = emergent fermions have no such direct connection to eigenvalues of Z. 



5 Free particle structures and counting on A4n 

In this section we extend aspects of our N = 2 discussion of the algebra of gauge invariant 
functions and the rings of scale invariant and gauge invariant differential operators to the 
case of general N. 

Following our considerations for the k = sector from Section 4, we show that the 
momenta of the free fermions are determined in terms of differential operators on A4.N- 
We show that the correspondence between states of the k = sector and the unitary 
matrix model (and hence free fermion wavefunctions) extends beyond two-point functions 
to a class of three-point functions. 

At N = 2 we have shown that the number of generators of the ring of gauge invariant 
functions on jMjy will be N 2 + 1. At N = 3, we will show that there is an interesting 
twist but that the above statement remains true in a refined form. 

Finally we study the m = n = k sector. This is the maximum possible value of k, in 
contrast to our studies of k = which is the minimum possible value. This sector consists 
of traces and multi-traces of Z^Z and we show that it may be mapped to N free fermions 
in a one-dimensional harmonic oscillator potential. This is a second, distinct appearance 
of free particles in complex matrix models. 



31 



5.1 The k = sector revisited 



We first observe that our construction of free fermion momenta as functions of differential 
operators in Zi,tij may be extended to general N in a slightly weaker form as follows. 

The construction in the previous section may be carried out for general N by identi- 
fying differential operators which measure higher order Casimirs. These will be traces of 
higher powers of the G\. We have not found closed form expressions analogous to (4.57) 
for higher N since this would require us to solve arbitrary order polynomials. However, 
since the Pi are integer or half-integer, they may always be determined in terms of the 
eigenvalues of the Hamiltonians [46], and hence implicitly in terms of differential operators 
in Zi,tij. We have thus identified an implicit map from k = operators to fermions on a 
circle for all finite N. 

We next conjecture that the k = sector may be descibed as the kernel of the differ- 
ential operator tr G 2 G 3 . Let us recall from Section 4.1 that the differential operator 



tr(G 2 + G 3 ) 2 = tr(G 2 2 + 2G 2 G 3 + G\) 



(5.1 



measures 6*2(7), an d so ti G 2 G 3 measures 



\ (C 3 ( 7 ) - C 2 (a) - C 2 (/3)) . 



(5.2 



Since for a k = operator 7 = (0, a, /3), we have that 

CM = C 2 (a) + C 2 (P) 



(5.3 



and so 




Since 



(G 2 y q Z) = 5\Z P 3 and - (G 3 )%Z^ = (5.6 



we have 



- trG 2 G3 {Z)^) = b\{Z^Z) k - (5.7; 



and since tr G2G3 acts via the Leibniz rule, the action of — trG^Ga on 




is that of the sum over all individual contractions 



n 




(5.9 



r=l s=l 
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Similarly the action of the laplacian 

1 = tr I 

. dZ azt 
on Z t ^ k l is given by 

n(Z)Z^)=5\5 k r (5.11) 

which is a Wick contraction using the two point function (1.1), and as before extends 
via the Leibniz rule. It was noted in [12] that the k = operators have no self Wick 
contractions and so we have 

nO k=0 (Z,Zi) = , (5.12) 

a result we shall use later in Section 6. For now we simply note that it is possible to 
construct simple examples which show that the k = operators do not comprise the full 
kernel of □. 

We expect however that the converse of (5.4) is true for any N 

tr(G 2 G 3 ) O = O = O k=0 (5.13) 

meaning that the kernel of tr(G 2 G 3 ) is exactly the k = sector. As a differential operator, 
tr(G 2 G 3 ) can be viewed as a modification of the laplacian which is invariant under scalings 
of Z and Z^. 

It is instructive to try and construct a counterexample to (5.13). From (5.2) we know 
that tr(G 2 G 3 ) O = is equivalent to 

C 2 { 1 ) = C 2 {a) + C 2 {(3). (5.14) 

The operator with labels 

a =[1,1], /3=[1,1], 7 = (fc = l j7+ = [l] >7 _ = [l]) (5.15) 

has Casimirs 

C 2 (a)=2, C 2 (/3) = 2, C 2 ( 7 )=4 (5.16) 

however this operator in fact does not exist since it fails our N = 2 constraint (4.16) in 
the form: 

d(a) + ci(/3) < N + k. (5.17) 

This example supports (5.13) and shows that it is sensitive to finite N constraints of the 
Brauer basis. 
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5.2 Three-point functions of k = operators 

In Section 4.6 we reviewed the map between k = operators and Unitary matrix model 
operators, and the result that the two point functions on both sides of the correspondence 
agree. Here we show that the same is true for a class of three-point functions. The 
correlators we consider are a subclass of the correlators 

(0 Al (Z, tf) A2 (Z, tf) 0\ 3 (Z, tf)) (5.18) 

where A\ = Pi Si is a short notation for the labels of the operators in the k = sector, 
and similarly for A 2 and A 3 . In terms of Brauer projectors the operators are defined by 

O r§ (Z, tf) = tx m , n {P RS Z <g> tf) (5.19) 

where P R § is defined in Appendix A. The subclass we study is that in which mi+m2 = m 3 
and ni + ri2 = 713. Performing the Wick contractions, we get 

(0 Al {Z, Zt) A2 (Z, Zt) As (Z, Zt)) = m 3 \n 3 \ tr m3jn3 ((P Al o P A2 )P A ). (5.20) 

The calculation is very similar to those in [2, 47]. It is convenient to express projectors 
as an integral over the U (N) group as 

P 1 = Dirwy J dU Xl {U ] )U (5.21) 

where Dirwy is the dimension of the U (N) representation 7; this follows from Schur-Weyl 
duality. We can therefore calculate (5.20) via 

tr m , n ((P^ o P A2 )P M ) = DimA 3 J dU 3XAa (Ui)ti ma , ns ((P Al o P M )U 3 ) 

= DimA 3 d Al d M \ dU 3 x Al {Uz)XA 2 (U 3 )xA z (Ul) 



= DimA 3 d Al d A2 g(A u A 2 ; A 3 ) 

= DimA 3 d Rl d R2 d Sl ds 2 g(Ri,R2;R3)g(Si,S 2 ;S 3 ) (5.22) 
where d A = d R § = d R ds and the following has been used to get the second equality: 

tr m , re ((P Al oP^ 2 )[/ 3 ) = d Al d A2XAl (U3)XA 2 (U 3 ) (5.23) 

and 

g(A 1 ,A 2 ;A 3 ) = [ dU 3XAl (U 3 )XAM)XA 3 (Ul) (5.24) 



is the Littlewood-Richardson coefficient which counts the number of A 3 in the tensor 
product A\ (g) A 2 . This derivation shows that integration over M. N can be done using 
Brauer algebras. 
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5.3 Finite N counting of single traces and multi-traces 

Since Z and do not in general commute, enumerating multi-trace operators in one 
complex matrix is an equivalent problem to enumerating multi-trace operators in two 
hermitian matrix models. 

Therefore an important check on the Brauer basis is that the counting of Brauer basis 
operators agrees with that of the counting of operators in two-matrix models [48, 49, 50, 
51, 52]. 

At finite N, we denote the counting of multi-trace operators built from two matrices 
with m of one type and n of the other as Q^ t (m, n). The following two expressions 
determine Q^ lt {m,n) in terms of group theoretical quantities [53, 48, 49]: 

Q% t (m,n) = Yl C{R,R,A)g{[m],[n]lV (5-25) 

R\-m+n,A\—m+n 
ci(R)<N,a (A)<2 

and 

(&(m,n) = E E 9{RuR2\ Rf- (5.26) 

Khm+n i?ihm i^hn 
Cl (R)<N ci(fli)<7V Cl (R 2 )<N 

Here C (R, R; A) is the multiplicity of the irreducible representation A of S m+n appearing 
in the tensor product of irreducible representations R (g) R of S m+n , and g(- , • ; ■) is a 
Littlewood- Richardson coefficient. 

Defining the finite N multi-trace generating function 

= 520%t(m,n) x m y\ (5.27) 

m,n 

the relation between the counting of single-traces and multi-traces is given by the Plethys- 
tic Logarithm [44, 45]: 

WW) = E^logK(^2/ fc ))- (5-28) 

k=l 

At N = 2 this leads to the five single traces identified in equation (3.37). 

At N = 3 the Plethystic Logarithm gives the single-trace generating function 

y^Q^ =3 (m,n) x m y n = 1 + x + y + x 2 + y 2 + xy + x 3 + y 3 

ra,n 

+ x 2 y + xy 2 + x 2 y 2 + x 3 y 3 - x Q y Q (5.29) 

The interpretation of this generating function is that there are 11 independent single trace 
operators along with one syzygy (algebraic relation) which occurs at order (m, n) = (6, 6). 
It would be interesting to find this relation explicitly. 



35 



We expect that the reduction multiplicities for B^ijn^n) to S m x S n , which we have 
denoted by M^p , satisfy the counting 

E( M ST = Q™t(™,n). (5.30) 

We show in the next section that this is true for N > m + n and it is also consistent with 
our calculations at N = 2 in Section 4.2. It would be interesting to prove it for general 
N, with or without using the connections to matrix models. 

5.4 Large N Brauer basis counting 

For N sufficiently large, i.eJV>m + n, we denote the counting of multi-trace operators 
by Q m t(m,n). The formulae in (5.25), (5.26) give rise to this quantity Q mt (m,n) when 
N > m + n. By Polya counting Q mt (m,n) is also given by (see [48, 49] & refs within) 

oo . oo 

n i _ (x r z yr) = E <Wm, n > m y n - ( 5 - 31 ) 

r=l ^ y ' m,n=0 

In equation (177) of [49] the following expression was derived: 

ci(l):EiCia)i=m c ; (2):E ! c i (2)«=n I i( J ' i( J ' 

The counting of Brauer basis operators is denoted N sb (m, n) and was shown in [12] to be 
given by 

N sb (m,n) = ^(M]) 2 . (5.33) 

In [12] it was argued that this formula correctly counts multi-traces at large N, namely 
that 

N sb (m,n) = Q mt (m,n). (5.34) 

In Appendices C and D we give two proofs of this fact, firstly by direct comparison to (5.32) 
and secondly by enumerating invariants in the reduction GL(N) x GL(N) — > GL(N). 

5.5 The m = n = k sector: Operators and free fermions 

We recall from the discussion at the start of Section 4 that the integer k is directly related 
to the minimum number of Brauer contractions involved in the terms which are summed 
to make up an operator in the Brauer basis. 
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For m — n — k, all terms in an operator involve the maximum number of contractions, 
which translates into the fact that these operators are multi-traces of the matrix Y = Z^Z. 
Since Y is hermitian we find the N fermions of the hermitian matrix model emerging in 
this sector, as follows. 

In this sector we have 7 = (k = m, 7+ = 0, 7_ = 0) and a = f3, so the projectors 
(defined in Appendix A) are in this sector labelled by a alone. We write 

p k=m = qj^ with 7 ag aboye> 

The projector is written in terms of the fc-contraction operator C( k ) defined by 

^( fc ) = zLj ^M 1 )! " '^(fc)fc' (5.36) 
aes k 

and the projector p a which projects the holomorphic half of V ® V to the represen- 
tation a. It is proved in Appendix E that the projector takes the form 

P « =m = C(fe)Pa (5.37) 

klDima 

and that the operator satisfies the following required properties: 

{p k= m) 2 = p k= m and ti k , k (P k=m ) = (d a ) 2 (5.38) 

where d a is the dimension of the S k representation a. The operators in the m = n = k 
sector therefore take the explicit form: 

ti Kk {P k=m Z m <g> Z*® k ) 

ti k , k {C {k)Pa Z® k <g> Z^ k ) 



klDima 



da £ ^AvCn ■ ■ ■ Cap^PaZ** ® Z*® k ) 



klDima 

= -^tr hik (C ir -.C k - kPa Z® k ®Z*® k ) 
Dima. 

= -J^Mp a Y® k ) (5.39) 
Dima 

where Y = Z^Z. So operators in the m = n = k sector are Schur polymonials constructed 
from Y. 

We may understand these results in the following way. First observe that Hl annihi- 
lates [Z^Zfji since H L = G 2 + G 3 generates the U(N) action on the lower index of 
and the upper index of Z, 

Z UZ, Z ] -> Z ] U ] (5.40) 
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and that the product (Z^Z)^ is invariant under this action. Traces of powers of Y are 
thus also invariant under (5.40). 

Hl measures 6*2(7) which implies that 6*2(7) = f° r an operators built from Y % y This 
is consistent with the fact that in the m = n = k sector 7 = (k = m, 7+ = 0, 7_ = 0) and 
so 6*2(7) = 0- We can consider a Casimir of the form tr(y^) 2 which measures the labels 
of the Young diagram. 

By the map discussed in Section 2.1, Schur polynomials in a hermitian matrix corre- 
spond to the states of N free fermions in a harmonic oscillator potential. The harmonic 
oscillator fermions observed here are a second emergence of free particles, distinct from 
those of the k — sector. 



6 Applications to integrable quantum mechanics 

6.1 Review of matrix harmonic oscillator quantum mechanics 

We now return to the dimensional reduction of M = 4 Super Yang-Mills on R x S 3 in the 
zero coupling limit, truncated to the sector of one complex matrix. As noted in Section 

3.2 one may choose the A = gauge while imposing Gauss's Law, yielding the quantum 
mechanics for the matrix Z(t) defined by the following action [54]: 

S = I ' dttT\ZZ ] - ZzA (6.1) 



It is well known that the holomorphic sector of the theory is equivalent to a system of 
non-interacting fermions in a one-dimensional harmonic oscillator potential [2, 8, 39]. As 
a subsector of M = 4 Super Yang-Mills extremal correlators in this sector are protected 
by supersymmetry [10, 11] and the states of this sector are dual to the LLM supergravity 
geometries [9]. 

Going beyond the holomorphic sector, we no longer have non-renomalization theorems 
so the connection to supergravity is not straightforward. Based on the following investi- 
gations, we will infer properties of any candidate string dual of the complex matrix model 
sector at zero coupling in Section 7. 

We first review the previous analysis of the above theory [2]. The momenta conjugate 
to Z l j and Z' l j are 

m = u Z i = z% ee n zt * = = z\. (6.2) 

Zj dZ) 1 z > dZt) 1 1 ; 

The equal time canonical commutation relations are 

[Z% n'J = i 8\5\ [Z\ nt'J = i 5\5\ (6.3) 
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so we can identify the conjugate momenta with matrix derivatives in the usual way using 
(2.5). We define the creation and annihilation operators: 



(6.4) 



Importantly, the dagger on does not signify hermitian conjugate of A. It signifies 
purely that this is a creation operator. The hermitian conjugate of A' % j is A J i . The 
canonical commutation relations become 

[A^A^}=S\5 k 3 [B),B^] = 5)5 k r (6.5) 

The Hamiltonian and U(l) current take the form 

H = trf-^^ + ZZM = tx{A^A + B^B) + N 2 
J = trTzA-ztAj = tr(A^A-B^B) (6.6) 

where N 2 is the zero point energy for N 2 harmonic oscillators in two dimensions. 

The ground state of this system satisfies A\0) = B\0) = 0. The corresponding (non- 
normalised) wavefunction \l/ = (Z,Z\0) is 

M Z i Zt ) = e~ tT ^\ (6.7) 

Holomorphic gauge invariant excitations of this system are defined by the constraint 
B\0) = and consist of operators built from A^ acting on the ground state. These may 
be written as 

ti n (a(AT m )\0) (6.8) 

where a is an element of S n , and controls how the indices are contracted to form either a 
single or multi-trace operator. A more convenient basis for operators of the form (6.8) is 
the Schur polynomial basis (for details see [2]): 

\®r) = Xr(^)\0) (6.9) 

where \R is the character of the U(N) representation R. Since 

we may write 

* R (Z,Zi) = XR (V2Z)e-^ zzt l (6.11) 
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This state has E = m + N 2 and J = m and is holomorphic in Z up to the exponential 
factor. If we triangularize Z and redefine the wavefunction by absorbing the Jacobian of 
the transformation into the definition of the wavefunction, it becomes a wavefuction for 
N fermions in the Lowest Landau Level of the Quantum Hall system [2, 8, 55]. 



6.2 Non-holomorphic sector 

The most general eigenstate can be constructed by acting with both A* and B' on the 
ground state, 

|*o) = O{A\B*)\0) (6.12) 

where 0(A\ B') is a gauge invariant polynomial constructed from m A^s and n B^s. 
The wavefunction of such a state may be written as 

Vo(Z t Z*) = (Z,Zt\* ) = 0(A\B^)e- tr ^\ (6.13) 

The Brauer Algebra may be used to organise the states above. Such states are analogous 
to those used in Section 4 and take the form 

\Kfw) = Ol^{A\B^) (6.14) 

where the labels are explained in Appendix A. This state has E = m + n + N 2 and 
J = m — n. 

Unlike for the holomorphic sector wavefunctions, we have 

0(A\ £t) e -tr(zxt) ^ {y/2Z, V2Zi) e ~ tTizzi) (6.15) 

because the derivative of Z inside A^ acts on Z which comes from the action of B^ on the 
exponential factor. For example we have 

ti(A^B^e~^ zz ^ = (2trZZt - N 2 ) e -*(^) (6.16) 

and in general the correct relation is 



$o(Z,Zt) = 0(A\B^)e-^ zz ^ = \e-?0(V2Z, y/2Z^] 



■tr(ZZt) 



(6.17) 



where □ is the laplacian tr ^g§f and the brackets indicate that the derivatives in □ act 
only on 0(y/2Z, V2Z^) and not on the exponential, e - ^ is defined by its series expansion; 
it was observed in (5.11) that the laplacian generates Wick contractions and so here e~ » 
performs a normal ordering, subtracting terms in which pairs of \/2Z and \J2Z^ have 
been contracted (c.f. [56]). 
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Note however that in a k = operator we have from (5.12) that 

□ O h=0 = 



(6.18) 



and so we can replace and with \/2Z and \/2Z^ respectively without worrying 
about the above subtlety. 

We can define operators corresponding to the G\ in (3.12) as follows. 



(Gi)) = (&B)) 
(G 3 )) = 



B^B\ 



(G 2 
(G 4 



- A^\A\ 



(6.19) 



Defining \A 1 -) = A l A0) and so on, using the commutation relations we find 

(G^IB^) = ^,|5ty (G 2 ))\A%) = P.\A*\) 

(G 3 ))\B%) = -5\\B^) {G 4 ))\A\) = -5\\A^) 



(6.20) 



which is the same as the adjoint action of the operators Gi defined in (3.12) on the 
matrices Z,Z^ (see equation (11) of [25]). 

The result is that we can define harmonic oscillator Casimir operators 



Ha= \Hi, H 2 , Hi, H 2 



Hr 



(6.21) 



by replacing Gj in (4.7) withGj. The eigenvalues of hatted Casimirs acting on Oa8-ij(-^> ^) |0) 
are the same as those of the corresponding unhatted Casimirs acting on ^ -{Z, Z^). 
This is because the same commutator manipulations can be done to evaluate both, and 
the arguments which prove that O2,0. i j(Z, Z') are eigenstates of the Casimirs in (4.7) also 
prove that ^ AA\ B^)\Q) are eigenstates of the hatted versions. 
We can take this one step further. Noting that 



zy e -S 



d 

d 
dZ* 



_ □ 
e 2 



(6.22) 
(6.23) 



and similarly 



[z%,«-!] 




( d V 


-\ 








\dZ ) 



e 2 



(6.24) 



then using (6.17) we derive 



A^^ (Z,Z^ = 0{A\ Bt) e" tr(zzt) 
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■tr(ZZt) 



■tr(ZZt) 



(6.25) 



where again the brackets indicate that the derivatives act only on 0(y2Z, y/2Z^) and 
not on the exponential. Similarly 



A)^o(Z,Z^ 



9 (£)',) ^ 



-tr(ZZt) 



(6.26) 



implying the following relation between G 2 and G 2 . 



(G 2 ))^ (Z,Z^ 



_□ 
e 2 



(G 2 ) 4 ? 0(V2Z, v^zJt 



-tr(ZZt) 



(6.27) 



Similar results apply to the remaining Gi, the Hamiltonians Hi as well as the canonical 
Hamiltonian 

H = H 1 + H 1 + N 2 = ti(A ] A + B ] B) + N 2 (6.28) 

whose action on wavefunctions ^(Z, Z^) can be written in terms of the (first-order) scaling 
operator H: 



H = Hi + H 1 + N 



N 2 . 



(6.29) 



Applying (6.27) and the corresponding relation for G3 we find that 



HV (Z,Z*) = e"T HO(V2Z, ^2Z^ 



■tr(ZZt) 



(6.30) 



A similar manipulation in the holomorphic sector was performed in Appendix A of [57]. 
Note that for a k = operator we have □ O k=0 = and so the above analysis gives 



H 



HO k=0 (V2Z, V2Z^ 



O k =\A\B^)e-^ zz ^ 
The inner product on wavefunctions may be derived using 

J [dZdZ ] ] \Z,Z ] ){Z,Z ] \ = 1 



-tr(ZZt) 



(6.31) 



(6.32) 



where [dZdZ^] = \\ i -dZ^dZ\-, as follows: 

1 



(*©il*0a) 



7T 



N 2 



[dZdZi] {O x (A\ B*)\Z, Z*)(Z, Z\0 2 {A\ £?+)> 



7T 



N' 1 



[dZdZ^] V 0l {Z, Z\)*Ch(Z, Z^) 



(6.33) 
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where ti n2 compensates for using non-normalised wavefunctions, and is found by imposing 

(*o|*o) = 1- (6-34) 
Using (6.17), the above expression (6.33) becomes 

(* 0l |* 02 ) = ^ j 1 {A\B^)e-^ zz '0 2 {A\B^)e- trZZ ' 

= -L J [dZdZ^] (e-SOi^,^)) (e-§0 2 (v / 2Z,v / 2Z t ))e" 2trZZt 

(6.35) 

and rescaling factors of two we have the result 

(*oJ*ft> = T^p / [ rf ^ f ] (e- n Oi(Z,Zt))( e - a 2 (Z,Zt))e- to ^ (6.36) 

which is the non-holomorphic generalisation of (A. 12) of [57]. 

In the next section we use the right hand side of the above equation to define an inner 
product on gauge invariant polynomials 0(Z, Z') rather than the harmonic oscillator 
wavefunctions ^o(Z,Z') which contain exponentials. 

6.3 Related integrable quantum mechanics models 

In the discussion above we related the action of the Hamiltonian H in terms of the (first- 
order) scaling operator H: We thus have an explicit map (6.30) relating the the action of 
H on its eigenstates 

Vo(Z,tf) = 0(A\B^e- tT{zzi) (6.37) 
to the action of H on its eigenstates 

0(Z,Z*). (6.38) 

The right hand side of (6.36) can be used to define an inner product on polynomial 
functions of Z, Z'. Explicitly this inner product is 



(0 1 (Z,Z^,0 2 (Z,Z^)) = J[dZdZ^}{e- a 1 (Z,Z^{e- a 2 (Z,Z^ 



e -tr zzy 



(6.39) 



Note that the two inner products ( • , • ) and ( • | • ) are defined on two different Hilbert 
spaces: 

( ■ , ■ ) : {Polynomials in Z, Z ] ) -> R 
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( ■ I ■ } : {Harmonic oscillator states} — > M. 



(6.40) 



By construction the inner products satisfy 

{0^,2^,0^,2^)) = (* 0l |* 02 ) (6.41) 
and each inner product is diagonalised by the corresponding Brauer basis 

Ol^jiZ, 2fi) and = Ol MJ {A\ Bt)|0). (6.42) 

The relevant Hamiltonians have the properties that 

• Ha are hermitian with the inner product ( ■ | ■ } 

• Ha are hermitian with the inner product ( • , • ) 

The inner product ( • , • ) can in fact be constructed by starting from the inner product 
arising from the zero- dimensional complex matrix model and requiring Hi to be hermitian. 
This involves subtracting Z, Z^ contractions and is discussed in Appendix F. 

In addition to the original Hamiltonian H, it is natural to consider the conserved 
charges e.g trG 2 , as Hamiltonians. Since these higher Hamiltonians were constructed to 
be simultaneously diagonalised with H in the Brauer algebra basis, we know these are 
solvable Hamiltonians related to the Brauer algebra. Related to these higher Hamiltonians 
are simpler ones which are obtained by replacing that hatted Cs with unhatted ones. For 
example H 2 = trG 2 . or H L = tr(G 2 + G 3 ) 2 define solvable quantum mechanics models. 
They are second order in derivatives as opposed to fourth order like H2, Hl, and hermitian 
in the inner product (6.39). 

Considering the expressions in (4.41) and comparing with equation (2.4) of [58] we 
see that these integrable quantum mechanics models are non-holomorphic generalizations 
of the Calogero-Sutherland model at a specific coupling (see also [59, 60] for related 
literature). A natural question is whether the Calogero-Sutherland model at generic 
coupling has such integrable non-holomorphic generalizations. Although we have not 
written out the Hamiltonians for general N explicitly as differential operators in terms of 
coordinates on A4. n it is clear that this can be done by changing variables from gl(N, C) 

to Zi, tij . 

7 Summary and outlook 

We described free particle structures hidden in matrix models of an N x N complex 
matrix Z. We related these structures to the geometry of the configuration space A4n 
of gauge- inequivalent configurations, a space of dimension iV 2 + 1. We showed that A4n 
supports an interesting class of functions, obtained from the gauge invariant functions of 
Z. The Schur decomposition gives coordinates useful for describing M. N . Integrals 
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over complex matrices give a measure of integration on A^jv which can be used to define 
an inner product on gauge invariant functions of 2^,t^. Following [12], Brauer algebras 
B^(m,n) give orthogonal bases which diagonalise the inner products which arise. Higher 
Casimirs constructed from the Brauer algebra, which resolve these orthogonal bases [25], 
give rise to a complete set of scale and gauge invariant differential operators on Ai N . 

Among the labels of the Brauer basis is a non- negative integer k. For any N the 
k = sector has states in one-to-one correspondence with those of N free fermions on 
a circle. These states correspond to the composite representations RS which play a role 
in two dimensional Yang Mills. The differential operators which measure Casimirs are 
polynomials in the free fermion momenta; for N = 2 we inverted these relations to write 
the momenta as algebraic functions of the differential operators. We conjectured that 
the k = sector is the kernel of a scale invariant version of the laplacian, the operator 
t^G^Gy. We also gave an equality of correlators between the unitary matrix model and 
the complex matrix model for a class of three-point functions of k = operators. It is 
important to note that while the usual emergence of free fermions in matrix models can 
be seen from a change of variables to eigenvalues [38], here the k = sector has no direct 
relation to the eigenvalues of Z. Indeed the operator t^G^Gy characterizing it involves 
derivatives with respect to Z% clS well as the off-diagonal Uj. 

Another interesting sector where states are counted by Young diagrams is the sector 
m = n = k. This is a sector of gauge invariant functions of (Z'Zy, which is the kernel of 
another second order operator on A4n, namely trG| as defined above (4.7). We observe 
that k appears to interpolate between radial and angular free particle systems on a plane. 
It would be interesting to further elucidate this in a stringy context. 

A precise understanding of the commutative ring of scale and gauge invariant differ- 
ential operators led us to computational results on the reduction multiplicities of repre- 
sentations of B]si{m,n) to S m x S n for N = 2. 

The connection between M. ^ and the space of gauge-invariant poynomial functions of 
is a more intricate version of the connection between WL N /Sn and symmetric polyno- 
mials. Likewise the connection between Brauer algebras and gauge-invariant differential 
operators on Ai n is a generalization of the connection between symmetric groups and 
differential operators on M. n /Sn- 

We present some avenues for future research : 

1. We wonder if some of these free particle structures can be obtained from the 
dual supergravity side of AdS/CFT in the sector which is SO (4) x 50(4) invari- 
ant. This would be a non-supersymmetric generalization of the LLM [9] discovery 
of supergravity geometries corresponding to the free fermions of the holomorphic 
sector of the complex matrix model [2, 8]. 

2. We have given explicit expressions for the free fermion momenta for the k = sector 
of the N = 2 matrix theories in terms of the original matrix variables. It is an open 
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problem to find explicit expressions for the coordinates of the fermions, and the 
wavefunctions as Slater determinants. It is also interesting to explore whether this 
would be useful for the computation of correlators. 

3. We have found a non-holomorphic generalization of the Calogero-Sutherland 
Hamiltonian at a fixed coupling. What is the physics of these non-holomorphic 
models? Can we observe Brauer Algebra wavefunctions in the laboratory? 

4. We have presented results on finite N counting of complex matrix model states in 
terms of Brauer algebras at iV = 2. These are related to reduction multiplicities 
for i?7v =2 (m, n) irreps into S m x S n irreps. What are these finite iV reduction 
multiplicities for general m, n, N, in particular for iV < m + n ? 

5. Our analysis has developed integrable quantum mechanics models for the 
space A4.N and exploited (Brauer) algebras to identify and organise interesting 
spaces of functions and differential operators on these spaces. A4n is a fibration 
over M, N /Sn which arises in hermitian matrix models and, like M. n /Sn, has differ- 
ent strata where the orbits qualitatively change their structure. While symmetric 
groups S n or their inductive limit Soo organise functions and differential operators on 
the symmetric product, the Brauer algebras BN{m,n) or similarly their inductive 
limit Bn(oo,oo) organise Mn- Results in matrix models, especially multimatrix 
models [61, 12, 49, 50, 51, 25, 62, 63] can give analogous results for other strat- 
ified spaces which arise as the space of inequivalent configurations. Is it possible 
to understand the role of algebras, integrable structures and hidden free particle 
systems intrinsically from the stratified geometries? What is the intrinsic charac- 
terization of stratified geometries which allow such structures? Studies of Hilbert 
space structures which mirror the strata in certain stratified spaces have been done 
[64]. Studying AAn from a similar point of view and finding its relations to the 
Brauer algebra description of functions and differential operators would be very 
interesting. 

6. There is a substantial literature discussing consistent truncations of the Maldacena 
duality. For example, it is known that the SU(2) sector defines a consistent trun- 
cation to all orders in perturbation theory [65]. Sectors such as the Z, sector 
are well-defined truncations at zero coupling. Assuming the strong finite N form 
of the Maladacena conjecture, and making the plausible assumption that consistent 
quantum truncations of a quantum field theory with a string dual have a string 
dual, we are led to ask: What is the gauge-string theory dual of one free 
complex matrix in four dimensions? Similarly what is the dual of the quantum 
mechanics from reduction of the Z, Z^ sector on M x S 3 7 For the large N Gaussian 
Hermitian matrix model (without the quadratic potential) there is the non-critical 
string considered in the old matrix model literature [27]. For double scaling limits 
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of the complex matrix models, we have the Type-0 string backgrounds [33]. For 
the large N hermitian matrix oscillator quantum mechanics, which is also a consis- 
tent truncation of the s-wave sector of M = 4 SYM in radial quantization, there 
is the proposal [66]. A well-known example of a duality between matrix quantum 
mechanics and M-theory is given by [67]. 

We do not have a clear answer to the last question, but the following remarks are 
suggested by the investigations in this paper. We conjecture that there exists a string 
dual of the matrix harmonic oscillator quantum mechanics discussed in Section 6 which 
has a 2 + 1 dimensional space-time and whose physics involves interacting strings and 
branes. The Zi coordinates are positions of N branes in 2 space dimensions. By analogy 
to the treatment in [68] we expect the variables Uj of the Schur decomposition to describe 
strings connecting brane i to j; here the triangular constraint (tij = for i > j) will make 
the dual qualitatively different from the standard system of strings and branes at weak 
coupling. This ought to be explained by an explicit construction of the string theory. The 
Hamiltonian H contains terms t-j^ along with Zi^-. Excitations involving polynomials 
in Zi have energies comparable to excitations involving t. This means that strings and 
branes have comparable masses. Usually string states have masses of order 1 (with l s = 1) 
whereas branes have masses of order l/g s . In this sense, the model at hand appears to 
have g s ~ 1. An interesting problem is to construct this strings and branes model in 
detail and to provide a physical interpretation for the labels of the Brauer algebra basis, 
in particular k, and their constraints at finite N. 
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A The Brauer algebra basis 



In this appendix we briefly introduce the Brauer algebra basis for gauge invariant poly- 
nomials in Z,Z\ The Brauer algebra B^{m,n) is used to construct a basis of these 
polynomials of degree m in Z and degree n in Z\ We will not need a precise definition 
of these algebras here, rather we will recall below how their representation theoretic data 
are used to label a useful basis. The definition of these Brauer algebras and their use in 
constructing a basis of gauge invariant operators is found in the original paper [12]. A 
more detailed review of the of the construction may be found in Section 2 of [63]. 

A Brauer basis operator is a linear combination of multi-trace operators built from m 
Z's and n Z^s and is written as 

Ol^{Z^). (A.l) 

These operators are constructed by viewing Z® m ® {Z*)® n as operators on V® m <8> V® n , 
composing them with elements in the Brauer algebra Q^p.ij and taking a trace [12]: 

The same construction can be done with the creation operators of the matrix quantum 
mechanics by replacing Z with and Z* with {B^) T where T denotes matrix transpose. 
These operators diagonalize the two-point function for Z, Z* at zero Yang- Mills coupling 
or the Fock space inner product for the states created by the A', of Section 6. 
The labels on the operator are as follows: 

1. a,/3 are Young diagrams with m and n boxes respectively, with c\(a) < A and 
ci (/?) < A. They label representations of U (A) as well as S m and S n respectively. 

2. 7 = (A;,7+,7_) where 

(a) k is an integer in the range < k < min(m, n) 

(b) 7+,7_ are Young diagrams with m — k and n — k boxes respectively, with 
Cl (7 + ) + Cl (7_) < A. 

7 labels a representation of the (walled) Brauer algebra Bpf(m,n). 

3. i,j are indices which run from 1 to the multiplicity M7g of the representation (a, 0) 
of C[S m x S n ] in the representation 7 of the Brauer algebra. 

The Brauer representation labelled by 7 = (A;, 7+, 7_) has an associated U(N) com- 
posite representation labelled by 7 C which is defined as follows. Using the usual notation 
in which a Young diagram with row lengths is written [7*1, r 2 , . . . , r^r], let 

7+ = [7*1, r 2 , • • • ,r p ], 7_ = [s u s 2 , ■ ■ ■ , s q ] (A.3) 
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then providing p + q < N, 7 C is given by 



7c = [ri,r 2 ,...,r p ,0,0, 




(A.4) 



where there are N —(p + q) zeroes inserted. In the language of the mathematics literature 
7 C is an iV-staircase with positive part 7+ and negative part 7_ [13, 15]. 

When we discuss Casimir operators we use the shorthand 6*2(7) f° r the ^ (N) quadratic 
Casimir of the representation 7 C , and similarly ^(7) or for the p-th row of j c . 

When k = the z, j labels are trivial and we have a = 7+, = 7_. Thus 7 is given by 
(/c = 0, a, /3) and so the k = operators are thus determined by two Young diagrams a 
and /3. To connect with the notation of 'coupled representations' in the two-dimensional 
Yang-Mills literature (see e.g. [69]), we rename 7 + — R,j- — S. If we substitute a unitary 
matrix in place of Z, the k = polynomials coincide with the 'coupled characters' Xrs{U) 
studied in the context of the string theory two-dimensional Yang-Mills [69] . 

At (m,n) = (1, 1), suppressing non-essential labels, the Brauer basis is 



Here we have suppressed 7+ and 7_ since for a k = operator it is always the case that 
a = 7 + and — 7_, and since for the above k = 1 operator, 7+ and 7_ are both the 
empty diagram. The multiplicity indices i, j are not relevant for this example. For further 
examples of Brauer basis operators see Appendix A.4 of [12]. 

B List of 7+ and 7_ at TV = 2 for given (m, n) 

Given (m, n) the possible 7+ and 7_ are listed below, along with 7 C as defined in equa- 
tion (A.4). Note that rj (and hence pj) distinguishes operators, as does r\. We use the 
shorthand (^2(7) for the U(N) quadratic Casimir of the representation labelled by j c . 



0^ m (Z,2fi) = trZtrZt-ltrZZt 
O^Z,^) = ^trZZt 



(A.5) 



(A.6) 
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List of 7+ and 7_ when m > n using d = m — n 



7+ 


7- 


7c 


k 


C 2 (7) 


[m] 


[n] 




[m, —raj 







m(m + 1) + n(n + 1) 


Iff) — i 


\n — 11 


[m 


- l,-(ra- 


1)] 


1 


(77?, — 1 I (77?, I -I- I 71 — 1 I ( Tl I 


[d + 1] 


[1] 




[d + 1,-1] 




n — 1 


(d+l)(d + 2) + 2 


[d] 







[d,0] 




n 


d(d + l) 


[d-1,1] 







[d-1,1] 




n 


(d-l)(d) 


Vdl \d\ 
1 2 1 ' L 2 J 







Vd~l \d\ 
1 2 1 ' L 2 J 




n 


ri i ( rn + 1) + Lfj ( lij - 



List of 7 + and 7_ when m < n using d = n — m 



7+ 


7- 


7c 


fc 


c 2 ( T ) 


[m] 


[n] 


[m, —n] 





m(m + 1) + n(n 


+ 1) 


[m — 1] 


[n-1] 


[m — 1, — (n — 1)] 


1 


(m — l)(m) + (n - 


-l)(n) 


[1] 


[d + 1] 


[1, -(<*+!)] 


m — 1 


(d+l)(d + 2) 


f 2 





[d\ 


[0, -d] 


m 


d(d + l) 







[d-1,1] 


[-l,-(d-l)] 


m 


(d-l)(d) 







d] d 
2 ' [2 


J [~<f 
2j ' 2 


m 


\l] ( |"f] + l) + [f 
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C Brauer counting at large N from Clebsch counting 

In this section we show that N s b(m,n) as defined in (5.33), 

N sb (m,n) = ( M a) 2 i (CI) 

j,A 



agrees with equation (5.32), 

TT — 

q(1)!q(2)! 



n < \ V tt (q(1) + q(2))! 

g m4 (m,n)= ^ ^ 11 „,mur9V • (C ' 2) 



Ci(l):E!Cl(l)I=m c i (2):E i c i (2)i=n J 

We first expand 



7 A 
min(m,n) 



n(m,n) / \ z 

E E E EE E»(*w«Wwfl 

fc=0 7+ h(m-fc) 7_h(n-fc) ohm /3hn \ Shk / 

= EEE 7+;aM*, 7-; £) ) (^(^7+;«M^7-;/3) J • (C.3) 

fc 7+,7- a,/3 \<5hfc / \5'hfe / 

Here (7 is the Littlewood-Richardson coefficient which is defined by 

g(5,~f + ;a) = — — — — ^ ^ XsWx-r+MXafa ° £r 2 ). (C.4) 

To simplify the expresssion (C.3), we use the orthogonality of characters 

v — v TV. 

2_^Xr(v)Xr{t) = S TatTT —- = S Ta ,T T Sym(T a ) (C.5) 



where T a is the size of the conjugacy class which contains a, and Sym(T a ) represents the 
number of elements which commute with cr: 

Sym(T a ) = c 1 (a)\l Cl{a) c 2 {a)\2 c ^ ) ■■■c n (a)\n Cn{<T) 

= Y[ci(a)li«W (C.6) 

i 

where c,(cr) represents the number of an z-cycle in a. Since a is an element of S n , we have 

Er=i^(°") =n - 

Using (C.5), some factors in (C.3) can be rearranged as follows: 

g ( 6 > 7+; 7+; «) 

7+ a 
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( fci( m - k)\ ) XsWx-r+MXafa ° ^2)X<5'(7-l)X 7+ ( r 2)Xa(ri o r 2 ) 



7+ « 



0"liO"2,Tl,T2 



r 1 ^ 


\k\(m - 


-k)\J 


ml 




(kl) 2 (m 


-k)\ 


m\ 




{k\) 2 (m 


-*)!. 



(m — k)\ 



XI ( r i ) xv in ) p^r i T 1 | \T T1 1 1 T T; 



m! 



CljC2,Tl,T2 



T1,T2 



p T2 | |T Tl or 2 
| X Xg(7"l)Xy(TL) | T \ \ T r 



CT1OCT2 I 



(C.7) 



Then we get the following expression for (C.3) 



min(m,n) 

N sb (m,n) = X 



m\n\ 



k=0 



(/c!) 4 (m - k)\(n - k)\ 



5,5' n,r 2 l J nor 2 | „, „„ |J, 



/'i I 



EL) ( fc ) E E 

fe V / V / ri,r 2 pi,p 2 



P1,P2 



Pl°P2 I 



E 



m\ n 



l 


l^nl 


I^pJ 




|2 IT- 

I I ± n °t"2 


KPl°P2 


l 


|r n | 





k \k) ^ ^ ^ \T T , OT A\T, 
k x 7 v 7 r l6 s fe r2e5 m _ fc p2e5 n _ fe 1 TlOT21 1 1 



n°P2 I 



E(l)(f)£ E E l T nl 2 |r„l|T w |— !— — i 



Sym(T Tl0T2/ 



n°P2 1 
Sym(T riop2 ) 



Sym(T T1 )S ym (T T2 ) S ym(T Tl )Sym(T P2 ) 

E E E En 



ciin) + ci(t 2 )\ ( ciin) + ci(p 2 ) 



The above can be rewritten as 

rain(q(l),q(2)) 

E E IIE E 

c l (l):J2lc l (l)=mc l (2):j:ic l (2)=n I k Ci (3):£ lc0)=k 



(q(1)-q(3))!q(3)!(q(2)-q(3))!q(3)!' 



We now compare to (5.32), which is the expression 



Q mt (m,n) 



E 



e n 



ci(iy.J2iCiWl=m Ci (2) : £ ;Ci (2)Z=n I 
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(q(1) + q(2))! 
q(1)!q(2)! ' 



(C.l 



(C.9) 



For any fixed cycle length in S m x S n consider the conjugacy class with q(1) cycles in 
S m and q(2) cycles in S n . The factor ~^r^?w^ is the number of ways of arrangements 
of (ci(l) + q(2)) objects with q(l) of one kind (say red) and q(2) of another kind (say 
blue). Suppose we lay our the objects in a line. We can take the first arrangemnet to 
be the one with q(1) reds on the left and q(2) blues on the right. Then we permute 
to generate the rest. A general arrangement will have q(3) blues on the left among 
q(1) — q(3) red objects and q(3) reds among q(2) — cj(3) blues on the right. Of these we 
have n? rrr^ r^vn x .., 2\ 2 ^' arrangements. Hence we get 

k\{ci{l)-ci(3)y. k\(c l (2)-c l (3)y. ° to 

(ci(1) + q(2))! q(1)! q(2)! 

q(1)!q(2)! q(3)!(q(1)-q(3))! q(3)!(q(2) - q(3))! ' 1 ' ] 

This proves the desired equality between (C.8) and (C.9), and so we conclude that the 
two countings (5.32) and (5.33) agree: 

N sb (m,n) = Q mt (m,n). (C.ll) 



D Brauer counting from GL(N) x GL(N) ->> GL(N) re- 
duction 

We now show that the Brauer basis correctly counts invariants under the adjoint U(N) 
action 

Z UZUl (D.l) 

We consider invariants under (D.l) constructed from objects of the form: 

Z\ ■■■Z im ; Z ]k \ ■■■Z^ kn , . (D.2) 

As far as counting invariants under U(N) action is concerned, the problem is equivalent 
to counting invariants under GL(N). 

The Lie algebra of GL(N) is just the full Matrix algebra M(N, C) and the symmetric 
algebra over S(M(N,C)) is decomposed into the direct sums [14] 

S(M(N, C)) = J2 V ^ ® (VwY ( D - 3 ) 

A 

as GL(N) ® GL(N) modules. The sum is over partitions with length at most N, i.e. 
Young diagrams with first column no longer than N. Restricting to S m (M(N, C)) leads 
to the restriction |A| = m, i.e we are looking at the case of Young diagrams with m boxes. 
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Decomposing the GL(N) x GL(N) into GL(N) we have 

S m (M(N, C)) = £ <?(t, 77; \)g(r, v- X)V V>V . (DA) 

T,T),V,X 

Here is a composite representation of GL(N). Therefore the set of invariants in 
S m (M(NX))®S n (M(NX))i is 

Inv I £ 9i,T,ri;\)g(T,u;\)g(T' ,ri';\')g(T' , 1/ ; \')V n>v ® V v >y > (D.5) 

v. r,r/,u,X X' ,t' ,T]' ,v' J 

which is nonempty only if 77 = 1/ , v — rj'. Hence the number of invariants is 

£ 9(r, V, %(r, u, X)g{r', u, \')g(r', r,, A'). (D.6) 

T,?7,i/,A,T',A' 



We relabel 



A — > a r — 7-7+ rj — > 5 

\'->P (D.7) 



to get 



a, /3, 7-l_ ,7- ,<5,<5' 



£ (X;s(7 +) 5;a)fl(7- 1 ^) J (£^(7 + ,5',«k(7-,5',/3)] (D.8) 

J,7+,7- \ S ) \ 5' J 



and using the definition of M2g (4.12) this becomes 



£(^) 2 (D.9) 



which is N s b(m,n) from (5.33). 



E Proofs for m = n = k projectors 

In this appendix, we shall show the operator (5.37) satisfies the following properties: 

( P k=m ) 2 = p k=m (R1) 

and 

ti k ,k{P k a =m ) = (d a ) 2 . (E.2) 
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E c<7(i)i • 




a£S k 




E ^ " " 


' C kk a 






E ■ ■ 


' C kk a 


a<=S k 





The second equation follows from the Schur-Weyl duality; 

7 

= ^ ^ yC(S fc xS fc ) ^ yB N (k,k)->C{S k xS k ) ^ ^ 

7, A 

In the second line, we have decomposed each irreducible representation 7 of the Brauer 
algebra into irreducible representations A of the group algebra of S m x S n . Acting with 
the projector p^ =m on this equation and taking a trace in V® k ® V® , we get 

tr k , k (Pt m ) = d(a,a) = (d a ) 2 (E.4) 

where we have used Divnq = 1 and M\ = 1 for 7 = (0, 0, k — m). 

The /c-contraction operator Ct k \ can be written in many ways, for example 

C( fe ) ■ 



(E.5) 



The second equality follows from 

aCq = C c{S) -jO (E.6) 
In order to show (E.l), we first calculate (C( k )) 2 : 

(C(fc)) 2 = E pC x \ ■ ■ ■ C h ip~ x aC x \ ■ ■ ■ Cfcfctr -1 

p,aeS k 

= E tT k(P~ 1(J )P C ll - --Ckk^ 1 

p,aeS k 

= Yl N c ^pC 11 ---C k - k a- 1 

p,udS k 

= J2 N c ^raC ll ---C k - k a- 1 

T,a£,S k 

= N k n k C {k) (E.7) 
where Q k is the Omega factor defined by 

Q k =J2 N c °~ k a (E.8) 
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where C a is the number of cycles in a. Using the equation (E.7), we can easily show that 
the projector (5.37) satisfies (E.l). 

We also have another interesting equation for Coa: 



which is a consequence of 



We finally prove (E.2) 



C( k )Pa = C( k )p a , (E.9) 



Cxi ■ ■ -C k - k a = Cii • • -C kk a ■ (E.10) 



trjb.fcCo-Cn ■ ■ -C k - k a 1 p a ) 



klDima 

d aGSk 

= ~~rt~^~ — tr fe,fc(Cll " " " C kk p a ) 

d' a 
Dima k a 

= — - — d a Dima 
Dima 

= K) 2 . (E.ll) 
F Constructing an inner product on polynomials 

The inner product on gauge invariant polynomials 0(Z, Z^) given in (6.39), 



(dCZ.Zt^Z.Zt)) = j^JidZdZ^fe^OxiZ^Z^fe^O^Z^Z^y 



(E.l) 



was introduced by identifying it with the integral representation of the inner product on 
matrix harmonic oscillator states \^/). 

In this appendix we show that this inner product can be derived by 

• Starting from the inner product arising from the two-point function of the zero- 
dimensional complex matrix model of Ginibre [1], 

{O x {Z,3fl),0 2 {Z,2fl)) a = j^J[dZdZ^]O^Z^)0 2 (Z,Z^e- trZZt 

(F.2) 

where the normalisation factor is the value of the integral with no insertions; 
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• Requiring Ha to be hermitian. 
We will see that this leads us to the inner product (F.l). 

The construction proceeds as follows. We know that H 1 , Hi, H 2 , H 2 , Hi have eigen- 
states given by the Brauer basis polynomials 0^q{Z,Z) with real eigenvalues. These 
eigenstates are a complete set of gauge invariant polynomials. So in fact any inner prod- 
uct diagonal in these labels 7, a, j3 

(0£a> OS*) = C^S aia2 S^ 2 (F.3) 

for some real f£ L a 1 wm guarantee that Hi, Hi, H 2 , H 2 , Hi are hermitian. 

We now give an explicit construction of such an inner product, which is also well- 
defined for general polynomials in Z, Z\ not just gauge- invariant ones. 

The basic idea is to define our inner product on degree 1 monomials in Z, Z^ and then 
extend to arbitrary monomials by Wick's theorem. 

{Z\,Z\) = 5 lk V l 
(Z%Z^) = 5% t 

(Z),Z^) = (F.4) 

We generalize to higher degree monomials 

(Z n - ■ ■ ■ Z im n Z^ v \ ■ ■ ■ Z tp " , Z kl , ■ ■ ■ Z km , Z^ r \ ■ ■ ■ Z tr " ) (F.5) 

by summing over different possible pairings of the Z on the left with the Z on the right, 
and the Z^ on the left with the Z^ on the right, with each individual pairing being given 
by (F.4). Very importantly we do not include contractions between pairs Z and Z both 
on the left or both on the right. In the Ginibre inner product (F.2), we have (1, ZZ} ^ 
which shows the inner product under construction is different to (F.2). 

To prove G 2 is hermitian with this inner product we first work with the basic pairing. 

(Z),(G 2 r q Z k l ) = b\{Zy,Z\) = 5 k q 5^5 3l 

((G a )*Z*,Z*) = V q {Z^,Z\) = 5^% (F.6) 

We thus find, on these degree 1 monomials 

i(G 2 y q ) h = (G 2 )%, (F.7) 

where h denotes hermitian conjugate. When we consider the action of {G 2 ) p q on a general 
pairing 

(Z h , ■ ■ ■ Z im - Z^ p \ ■ ■ ■ Z tp " , (G 2 ) p Z kl , ■ ■ ■ Z km , Z* r \ ■ ■ ■ Z tr \ ) (F.8) 
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we can use the fact that {G2) p q acts as a derivation, so that the right factor becomes a 
sum of terms with the {G2) p q acting on each successive Z ov Z'. The action on Z^ gives 
zero. For each term in this sum, the inner product is a sum over Wick contractions. For 
each Wick contraction of the form 

(Z,GZ)(Z,Z)---(Z\Z^)--- (F.9) 

we can move the {G 2 ) p over to the left to give (G 2 ) % using (F.6). We can recollect the 
sum over Wick contractions to get 

{{G 2 )%Z\ ■ ■ ■ Z^ Jm Z^ qi ■ ■ ■ Z^Z k \ ■ ■ ■ Z k - lm Z^ Si ■ ■ ■ Z^J. (F.10) 

This establishes for any monomial in Z, Z* that ({G 2 ) p q ) h = {G 2 ) q p , and by linearity this 
extends to any polynomial. Having established 

{{G 2 y q ) h = {G 2 y p (F.n) 

it easily follows that 

(trG 2 f = trG 2 and (tiG\) h = trGl (F.12) 
and similarly we find 

(((''■:)',)'■ = (G,) J , 

(tiG 3 ) h = trG 3 

(trGl) h = trGl 

(ti(G 2 G 3 )) h = tr(G 3 G 2 ) = tr(G 2 G 3 ) (F.13) 

where the last equality follows since the entries of G 2 and G 3 commute. 
We can also derive the above relations by noting that 

W " {ml- 

(^rS = Ui = Z\ {jz) r = {G 2 )' p (F.14) 

and similarly for G3 etc. 

The above proofs work by construction since we have defined our inner product to have 
the same properties as the oscillator inner product for A^B^ and exploited the similarities 

(G 2 )< ~ (A^A)), (G 3 )) ~ (B^B)). (F.15) 
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We now derive an integral form of this inner product. We begin with (F.2) and normal 
order by removing all contributions to the inner product from self-contractions in the 
wavefunctions. It was observed in (5.11) that the laplacian generates Wick contractions 
so we define (c.f. [56]) 

□ 2 

:0(Z,Z j ): = (l-D + — + ---)0(Z,Z^ = e- n O(Z,Z^) (F.16) 

and our inner product becomes the following modification of (F.2): 

(0 1 (Z,Z^),0 2 (Z,Z^)) = (:Q 1 (Z,tf):,: Q 2 (Z,tf):) G 

J [dZdZi] (e^ 0l (Z,Z^(e- n 2 (Z,Z^ 



(2vr) J 



e -trZZt 



(F.17) 



which is (F.l) as we set out to show. 
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